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ABSTRACT 

In its simplest form, a radial face mechanical seal 
consists of two annular rings, one of than is moxmted 
on a flexible support and is called the primary seal 
ring and the other is set on the rotating shaft and 
is termed as the secondary seal ring. For smooth 
operation of face seal over a long period of time, the 
faces of the seal have to remain properly separated* 

Thus an adequate amount of separating axial force is 
required to be maintained for this purpose. But too 
big a separation leads to massive leakage of the 
lubricant whereas too small a separation remits in 
the rapid wear of the faces . Thus short seal life and 
sudden seal failure, because of uncontrolled leakage 
or wear or both are common problems. 

Recently it was held that restoring moment (the 
moment about x axis) and transverse moment (the moment 
about y axis ) have direct bearing on the stability of 
face seals. Both hydrodynamic and hydrostatic effects 
were shown to have enormous influence on various seal 
characteristics namely, axial force, radial force, leakage^ 
frictional force, transverse moment and restoring moment. 
Misalignment and coning of the faces is shown to be inevit- 
able phenomenon and they influence these characteristics 
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considerably. Also, squeeze effects are shovm to have 
significant effects on various seal characteristics* 

The non-Newtonian character of the lubricant was 
shown to have insignificant influence on the performance 
of a radial face seal. But these findings were made 
without taking into account the inevitable misalignment 
of the faces and the flexibility of the face material. 
However, it is likely that the study of radial face seals 
V7hen coupled with the concepts of misalignment and flexibi- 
lity may lead to significant non-Newtonian effects. 

Further, it is likely that the study of these non -Newtonian 
effects through the power law model may yield additional 
information which will help in the designing of seals. 

The present thesis is designed to study some of 
these aspects of a radial face seal_. The seal faces 
are assumed to be narrow in that it renders Reynolds 
equation worth yielding analytical solution for some 
cases, as a consequence of elimination of the circumfer- 
ential pressure gradient frcro the original Reynolds equation. 
In fact, for narrow seals, the circtmferential pressure 
gradient is neglegible as compared to the radial pressure 
gradient. The present thesis is devided into two parts ~“ 
Newtonian effects and non-Newtonian effects. The 
chapterwise brief discussion of the work presented 


is as follows 



vi 


Chapter I is indroductory type. Starting from 
the basic definition of seals, this chapter deals with 
the various aspects of designing of seals, their 
applications and their classifications. Besides, the 
chronological davelofxnents in the field have also been 
touched upon. Some of the important investigations 
have been described in detail. 

Chapter II deals with the effects of centrifugal 
inertia on the performance of a parallel face seal. 

It is found that both axial force and side leakage increase 
with an increase in inertia parameter whereas transverse 
moment decreases. 

In Chapter III, squeeze effects on face seals are 
studied. The faces of seal are considered to be misaligned 
and one of them carries coning on it. It is concluded 
that squeeze with coning may be desirable where the large 
values of separating forces produce massive leakage, 
because coning causes reduction in these forces. 

In Chapter IV, hydrodynamic effects in misaligned 
radial face seals, with porous faces, are analysed. It 
is concluded that when there is a case of excessive leakage, 
due to large values of separating forces, rendering the 
faces porous reduces the leakage thus saving the face seal 
from early failure. 
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In Chapter V, the problem of misaligned radial 
face seal is dealt with where one of the faces is 
onsidered to be flexible which may be due to coating 
of a soft material on the face seal. It is found that 
the axial force increases with the elastic parameter 
whereas the radial force decreases. The conclusion is 
that the flexibility of faces improves the stability 
of face seal. 

In Chapter VI, the problem of misaligned radial 
face seal is considered, where the non-Newtonian behaviour 
of the lubricant is accounted for through the power law 
lubricant . It is seen that the axial force and the 
tilt moment are increased by the use of pseudoplastic 
fluids as Ixibricant. However, they are decreased by 
the use of dilatants. 

In Chapter VII, the effects of misalignment on 
radial and frictional forces are studied. It is 
assiamed that the lubricant obeys the power law. The 
curvature effects neglected by earlier workers are 
also given due attention. It is seen that the pseudoplastic 
fluid increases the radial and frictional force whereas 
the dilatant fluid decreases them. Furthermore, higher - 
the radii ratio higher the radial and frictional forces. 
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In Chapter VIII, elastohydrostatic Ixibrication# with 
power law lubricants, of circular plate thrust bearing 
of the type used in face seals is studied. It is found 
that an elastic layer on the surface improves amply 
the bearing perfomance. When extended to face seals, it 
amounts to saying that axial separating forces are increased 

ft 

by an elastic layer on the faces of a face seal, 

A part of the work embodied in the thesis has been 
published in the form of the follov/ing research papers, 

1. "Squeeze Effects in a Misaligned Radial Pace Seal 
With Coning", Wear, Volume 85, 1983, page 143, 
(Chapter III). 

2. "Radial and Frictional Forces in Misaligned Radial 
Face Seal With Non-Newtonian Fluid", Presented 

^ at The Winter Annual Meeting of the American Society 
of Mechanical Engineers, Phoenix, Arizona, November, 
1982, (To be published in Jour. Fluids, Engng. , Trans . 
ASME). (Chapter VII). 

3. ''Elastohydrostatic Liibrication of Circular Plate 
Thrust Bearing With Power Law L-ubricants " , Journal 
of Lxibrication Technology, Transaction of American 
Society of Mechanical Engineers, Volume 104, 1982, 
page 243, (Chapter VIII), 
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film thickness corresponding to the outer periphery 
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pressure in the porous shell 
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f 

measure of the f lexibility, i.e. s ss - 

E' 

velocity in r direction 
velocity in y direction 
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a measure of the porosity of seal faces, i. e.a=i2 ^ 
angle of coning 
angle of tilt 

non-dimensional coning parameter, i.e. 6 = 

7r ^ 

2 

non-dimensional tilt parameter i.e» e = 

inertia parameter depending upon angular velocity 

poisson's ratio 

shear stress 

permeability 
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CHAPTER I 

GENERAL INTRODUCTION 

In any lubricating system the purity and the 
maintenance of the quantity of lubricant is highly desirable. 
'Sealing in' the lubricant and 'sealing out' the contaminating 
foreign material are the measures that serve this purpose. 

The function of seal is to separate pressurized fluid where# 
for instance# a moving shaft passes through a machine 
housing or passes from one part of the machine to another. 

This requires prudent designing of the seal arrangements 
and the selection of the seal material. In the selection 
of any combination of seal arrangements and material the 
objective should be to obtain a design which will show the 
least amount of wear in service. It is also to be stressed 
upon that the designing may be of such a form that whatever 
wear takes place should be on the seal and not on the shaft. 

In other words the amount of wear on the lubricating system 
should lessen on account of sealing. Por example# a leather 
seal will wear depending on the quality of the material of 
seal. But the wear on the seal can be reduced in those 
design where it is possible to fit a spring with the seal 
mounted on it to hold the seal close to the shaft. 

1.1 SEAL DESIGNING 

Since the lubricant in a mechanical system may have 
to operate under certain conditions# the arrangements of any 
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seal for the mechanical system calls for the consideration 
of all such conditions by a designer. Some of the important 
conditions are discussed as follows- 

1*1«1 Temperature ! It is well known that higher temperature 
of the lubricant results in higher fluidity. Thus if the 
seal is not sufficiently tight the leakage of the lubricant 
will be enhanced. Thus a seal designed for a low-temperature 
operating condition will fail to sufficiently seal in the 
fluid when the lubricant has to fionction under high temperat- 
ure. The seal designer, therefore, considers the operating 
conditions when choosing the seal material and planning 
the spring arrangements, for high temperature may adequately 
influence the characters of the material used. 

1.1.2 Pressure ‘ Under most conditions of service in a 
pimp or gearhousing or a roller bearing some pressure will 
prevail within the lubrication region. This means that 
in operation there will be more of a tendency for the 
lubricant to leak out than for dust or dirt to be drawn 
in. On shutdown, however, a partial vacuum may develop 
in some lubricating system upon cooling and cause dust 
laden air to be drawn in. In view of this ,it is not alv/ays 
possible to seal all types of bearings outright against 
the leakage or the infiltration of dust or foreign matter. 
Thus some designers make additional arrangements for more 
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ef f ect«^ overall performance of the systsn. The use 
of external f lingers is one of such additional arrangements. 
External f lingers are useful where water may run down 
(or along) the shaft. They are also helpful in providing 
security to the bearing seal against the wear-borne effect 
of contaminating particles . 


1.1.3 Wear * To design a seal that lessens the wear it is 
required that the material comprising seal-faces should 

be properly chosen. In the case of unsuitable combination 
of face material the wear may reach disastrous proportions. 
In practice, a badly designed seal which has a good combin- 
ation of face material is often preferable to a better 
designed seal with unsuitable face material. The porosity 
of the material, for instance, is harmless as long as the 
pores are not interconnected. Such hollow pores have the 
advantage that they may store lubricants and may therefore 
have improved perfomiance. Examples of such materials 
are metal oxides, steel containing graphite, and castings. 

A similar effect is produced through the action of wear 
on unhomogeneous materials with particles of varying 
thickness in the matrix or differing in thermal conductivit- 
ies and expansion coefficient. 

1 . 1.4 Type of lubricant I The type of lubricant to be 
used sufficiently influences the kind of designing a seal 
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really requires. Fluid oils operating at temperatures 
within average limits (from 50° to 150°C} require tighter 
sealing than greases which are made from oils of approximately 
the same viscosity. The soap content in the average bearing 
grease has a certain amount of restraining effect in itself. 

Thus the designer of the seal must have beforehand 
sufficient information regarding the type of lubricant to 
be used, viscosity of the lubricant , the operating temperature, 
the dust density of the surrounding or the other contaminating 
outer materials, the operating speed of the system, the kind 
of products that might be ruined by the lubricant in case 
of its leakage, e.g. textiles by the lubricants, etc. It is 
also important to remember that for any seal to be effective 
the surfaces against which it rvibs must be smooth and well 
polished. 

1,2 TYPES AND CLASSIFICATION OP SEALS 

In Industry a large number of seals are used to meet 
the many and varying problems presented by fluid sealing 
needs. Depending on the operating conditions and requirements, 
gaskets , 0-rings and welds are used for sealing stationary 
machine joints J soft packings, moulded elastomeric rings, 
piston rings, metal bushings and diaphragm seals are used for 
backward and forward movements. 
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Seals mounted on springs, stuffing boxes, lip seals, 
bushings, spiral groove seals and mechanical seals made out 
of a very large number of materials are the systems that 
are commonly used for sealing rotating shafts. Bearing 
engineers, group seals into the following three broad 
classifications according to the principles of design 
involved!! 

(1) Those involving actual frictional or sliding contact 
with the shaft. These require a relatively flexible 
sealing material such as leather ,f elt, cork or some 
synthetic material which can be pressfitted into 
the closure so that light contact prevails with 

the shaft. 

(2 ) Those involving grooves or metallic devices such 

as stuffing boxes which require a certain amount of 
running clearance between the shaft, and make use 
of capillary attraction to keep a film of lubricant 
in the clearance space. 

( 3 ) Those which utilize the effect of centrifugal force 
by providing flingers vrhich will keep dirt from 
entering the housing. 

When the leakage flow is radial between plane 
surfaces the mechanical seal is termed as "radial face seal'" 
and when the leakage flow is axial between the cylindrical 
surfaces the mechanical seal is termed as 'bushing seal . 
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Axial forces control the leakage in radial face seals 
and the size of the radial clearance between bushing and 
shaft controls the leakage in bushing seals. The face 
seal is more and more replacing stuffing boxes for rotating 
shafts, having the advantages of smaller leakage losses, 
reduced maintenance and longer life. 

1.3 FIELD OF APPLICATION 

It was the automobile industry that recognized the 
advantages of mechanical seals and first used them in the 
water pumps of internal combustion engines . Today automobile 
water pumps are all fitted with face seals. After face 
seals had proved themselves in the motor industry, the 
pump manufacturers and the large chemical firms began to 
install an ever ■ -increasing ntamber of mechanical seals in 
a rapidly expanding range of applications. The petroleuiti 
industry went so far as to prescribe the exclusive use of 
mechanical seals in some refineries, pumping stations and 
tankers for reasons of safety and economy. 

In the motor industry, face seals are used in 
^ abundance. They serve to seal crankshafts, water pumps 
and auxiliary pumps in motor vehicles, diesel engines and 
ships -engines. They seal bearings in tracked and cross- 
country vehicles, preventing the extrusion of grease and 
oil and the ingress of foreign bodies or dirty water. The 
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operating requirements for motor vehicles lie mainly 
in the low and middle range of pressure, speed and load- 

Face seals are frequently used in household 
machinaries such as v/ashing machines, driers, spin-driers, 
dish-washers and meat-mincers, as well as pumps for oil 
heaters and hot water distribution. Power and pump 

e- 

industries are also greatly benifited by the use of seals. 
Here, seals are used in water turbines, boiler feed 
pumps, gas circulators and to an increasing extent in 
nuclear reactor technology on high pressure rotary pumps 
and other control rod mechanisms. Furthermore, fans, 
compressors, refrigerators and vacuum pumps are fitted 
with face seals. 

In the petroleum and chemical industries mechanical 
seals find wide applications and are required for the 
shafts of agitators, sprinkler systems and mixers, steam 
heated calenders, centrifuges, blowers and pimps for 
varied chemical products. In aircrafts and rockets, 
mechanical seals are used for sealing gas turbine shafts, 
turbos erper chargers, hydraulic units and booster pumps 
for fuel and liquid gas. Here the rubbing speed can 
be very high. Apart from this, seals find applications 
in many other industries. 
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1,4 EXPLORATIVE DEVSLOfMENTS 

The first simple face seal appeared in about 1900 
when very little information about functioning of seals 
was available. Although the use of seal dates back to a 
period far more prior to 1900, but this simple face seal 
was an improvement over all previous face seals e.g. trad- 
itional shaft packing and stuffing box. 

Over the past several years advances in contacting 
mechanical face seal technology have been promising, but 
considerable work still needs to be done. Further improve- 
ment in seal design to increase seal life and reliability 
would result in a sufficient cost saving to the users. 

Also, it is known that the amount of energy consumed for 
pimping purposes due to seal friction losses can be signi- 
ficant. With the present state of the art of contacting 
face seal design, it is difficult to design a seal with low 
leakage, low friction, high reliability and long life. 
Considerable explorative work has been done over the past 
three decades. Judging by the progress hitherto one may 
expect continued exploration and increasing knowledge 
leading to evolution of still better seal design technology. 

The earliest available literature dates back to 
1946 when Pogg^s experiments [ 1 ] on simple parallel 
face thrust bearings showed that considerable load could be 
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carried by flat faces without physical contact. Four 
years later, in 1950, Brkich L 2 ] found out that the 
surface tension acting at the interface between liquid 
and air, at the low pressure edge of the seal faces, 
created a pressure differential to contain the liquid. 

The sealing pressure so created can be computed knowing 
the surface tension and assuming that the maximum gap 
between the seal faces is equal to the sum of the waviness 
of two face surfaces. In 1950 itself Salama [ 3 ] investi- 
gated that surface roughness produced by machining has 
significant bearings on pressure producing mechanism. 

In 1954 a paper by Wood L 4 ] appeared in which ejqjloration 
about seal-life was discussed. In conservative refinery 
practice, typical seal-life was noted as being about 
2000-8000 hours. In domestic water pumps, life is probably 
of the order of 15000 hours, i.e. two years of continuous 
service. It was also investigated that for very high 
face velocities and poor lubrication, life may be only 
a matter of minutes. In 1957 Zienkiewic^ L 5 ] found 
that for some configurations, thermal and pressure defor- 
mations of the bearing surface produce a convergent film 
vrhich is load bearing. Surface tension as the means 
of sealing was discussed by dagger [ 6 ] in the same year. 
He showed that the maximtim pressure that could be sealed 
by this method was in the range of 2 5 to 50 psi for common 
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fluids assimiing a 20 microinch gap between the seal faces. 
Whitley [ 7 ], in 1959, used grooves and slots in the seal 
faces to generate fluid pressure to maintain a full fluid 
film and thereby supporting load and minimising friction. 

Thus, although 1950s witnessed a considerable work 
on face seals, yet the theory of face seal did not progress 
to the point of providing a general useable design of face 
seals. Still more investigations followed in 1960s especially 
after Denny L 8 ] had carried out some experiments with 
plain thrust surfaces of the type used in radial face seals. 

In the range of conditions covered by these experiments, 
fluid pressures between the rotating annular surfaces were 
always greater than those predicted on the basis of laminar 
flow between parallel surfaces. He inferred that the factors 
contributing to this wore as follows • 

(a) Fluid pressures were generated between the faces 
by the shaft rotation. 

(b) Wear of the faces tended to produce a convergent 
gap in the direction of flow. 

Both effects resulted in stabilising action to 
prevent the faces from separating inadvertently. It was 
tentatively suggested that the generated pressures were 
due to vibration of the faces and the wear was due to impur- 
ities in the liquid. A centripetal pumping action of 
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^PPrsciable magnitude was shown to occur under some 
conditions, so that the flow against a pressure gradient is 
possible . 

Nahaldi and Osterle [ 9 ] studied the effects of 
face misalignment using different analytical techniques 
and found that axial vibrations increased or decreased the 
load bearing ability depending upon particular example. 

In 1960 Tanner [ 10, ll]shov/ed that the fact that the 
average pressure over the seal faces was greater than the 
oil pressure as was pointed out by dagger [ 6 ]was unlikely 
to be explained by the effects of non-Newtonian behaviour 
of the fluid. Instead, as he suggested, it is possible 
that elastic deformation with this factor and surface 
tension may give a different result. He further mentioned 
that non-Newtonian effects do exist but the small order 
of magnitude of these effects will render them unimportant. 

In 1961 a ntamber of papers [ 12-17] were presented in BIIRA 
International Conference on Fluid Sealing. Both theoretical 
and experimental investigations were performed to have insight 
into the effects of frictional forces and surface roughness 
on sealing. Here the findings by Mayer [ 15 ] shares quite 
a significant contribution in that he found that higher 
pressure could be sealed with balanced seals than with 
unbalanced seals. This originally raised the question of 
why relatively small changes In balancing can alter the 
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maximum sealing pressure. This difference could be 
ascribed to the different lubricating regimes that might 
have existed between the seal faces. 

Whitelay [ 18 ], in 1962 , brought out that the 
performance of seals is improved by using grooves. In 
1963 Davies [ 19 ] and Lyman and Saibel [20 ] considered 
the effects of surface deformation and face misalignment 
respectively on leakage and axial forces. In 1967, 
some papers [ 2 1-2 3 ]_ presented in the Third International 
Conference on Fluid Sealing [ BHRA ] dealt amply with 
the hydrodynamic aspects of face seal lubrication. The 
pressure distribution in high pressure face seal was also 
' touched upon. In 1963, the papers bearing the reference 
[ 24-27 ] discussed the hydrostatic and hydrodynamic 
large scale waviness and microasperities. The papers 
[ 28-35 ] published in 1969 further considered misalign- 
ment and vibrations of seal faces and microasperities 
of seals and demonstrated that these can cause hydrodyna- 
mic load support. 

The period during 1970s can be cited as to have 
drawn appreciable amount of attention of various research- 
ers to the functioning of a mechanical face seal. The 
effects of waviness, coning, misalignment, flexibility, 
roughness on leakage, axial and radial forces and moments 
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were mainly studied. Gibson et.al. [ 36 ] studied, 
in 1970, the elastohydrodynamic problem of reciprocating 
rubber seals for hydraulic actuators. The problem was 
that of a lubricated highly elastic surfaces in pure sliding- 
In 1971, Iny et.al. [ 37,38,39 ] presented a comprehensive 
theory of sealing based on the principles of hydrodynamic 
lubrication with the faces separated by a fluid film. 

He established that the hydrodynamic theory can be applied 
to the design of mechanical seals for more severe operating 
conditions. Field and Nau [ 40 ] in 1972 measured exper- 
imentally the oil film thickness and pressure distribution 
under an actual rectangular piston seal. He concluded 
that there is considerable difference in the behaviour of 
0-ring seals and rectangular seals. Lohou and Godet [ 41 }, 
in 1973 found that misalignment in radial face seal can 
cause hydrodynamic fluid pressure load support. Stanghan- 
Batch et.al. [42 X applied hydrodynamic lubrication 
theory to an idealized v;avy face seal. Extending the 
theory to cover radial flows in the fluid film the author 
indicated that the hydrodynamic film pressure combines with 
the curvature of the seal faces to induce an inward flo^i/ 
of fluid across the faces. The frequently observed inward 
flow phenomenon may be utilized to assist sealing. Under 
certain circumstances, however, the inward pumping effect 
can be a cause of seal failure, particularly at low fluid 
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pressures, and the choice of an inside or outside seal 
arrangement should be made to suit the duty required. The 
results of his work suggest that mechanical seal faces 
should not be lapped flat, but one of the faces should be 
given a deliberately wavy surface. 

In all conditions of seal functioning the hydrodynamic 
effects are not the only ones that maintain lubricated gap 
in the face seals. This has led to the study of hydrostatic 
load support. Its significance impelled several researchers 
to use new techniques of predicting and controlling the 
radial face seal profile . Snapp and Sasdelli [ 43 1 assumed 
that the hydrostatic pressure drop across the seal faces 
is highly influenced by the gross deformation of the seal 
rings. In 1973 itself Metcalfe[ 44 ] carried out his vrork 
on pressure caused effects and found that pressure caused 
deformation resulted in divergent seal profile. An account 
of various developments have been embodied into the book 
by Mayer [45 ]• Radial variations in seal clearance affects 
both hydrodynamic and hydrostatic behaviour of primary seal 
ring. Moreover, the existence of full fluid film depends 
on both hydrodynamic and hydrostatic components of pressure. 
These investigations have been made by Zuk [ 46 3/ Shapiro 
and Colsher [ 47 ] and Haardt [ 48 1 in 1974. In the same 

Ruskell and Westcott [ 49 ] gave an accoxant of analysis 


year 
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made on the performance characteristics of thin lipped 
metal seals as compared to that of piston-ring seals . The 
work carried out in 1975 bears the reference L 50-53 ] 
which, apart from discussing the seal dynamics , involves 
the effects of vibrations of seal ring in one degree of 
freedom. Field and Nau [ 53 ] also presented solution 
of elastohydrodynamic equations based on iteration scheme 
for a rectangular rubber seal. 

In 1976 Ludwig [ 54 ] , Ludwig and Allen [55 ] 
and Metcalfe [ 56 ] studied further hydrodynamic pressure 
load support and pressure and thermal caused distortions. 

It was found that thermally caused rotation from friction 
results in a convergent seal profile and pressure caused 
deformation results in a divergent seal profile. Thermal 
deformation was further studied by Li [ 57 ] and Kilaparti 
[ 58 which led to the inference that thermal deformation 
prevents seal faces from remaining flat. To have better 
prediction of seal performance the consideration of face 
misalignment was deemed to be necessary. Ruskell [ 59 ] 
made a significant comparision between theoretical and 
e 3 <perimental results for elastohydrodynamic lubricauion 
of metal seals. It was demonstrated that an unusual seal 
profile is caused by the deformation of the seal. It was 
also predicted that thermal effects, non-Newtonian behaviour 
and thermoelastic effects have possible relevance to seals. 
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As in the ease of microasperities and misalignment/ 
large scale waviness can cause hydrodynamic fluid pressure 
load support. This was amply supported by the investigations 
of Lebeck et.al [6l]and Lebeck [ 61/62 ] in 1977. Lebcck 
et.al [ 63 ], in 1978, gave experimental account of the 
effects of asperity on leakage and seal profile. Etsion 
[ 64/65 ] showed that the effect of curvature of a radial 
face seal is negligible. The effect of non -axi symmetric 
hydrostatic pressure duo to radial variations in the film 
thickness is considerable, which in turn influences seal 
stability and wear of the outer and inner diameter. Hughes 
et.al [ 66 ] modeled a seal with a phase change using 
an idealized heat transfer model. Their results showed 
that the hydrostatic load support in a two phase seal is 
greater than for an all liquid or all gas seal . Patir 
and Cheng [ 67 ] suggested for incompressible flow that 
the f^ow equation should be multiplied by a factor Q 
depending on operational conditions. Metcalfe et.al [ 68 ] 
indicated that to improve seal performance in some cases 
the flatness of faces is deliberately disturbed. Wobble 
observation of primary seal ring is reported in [ 69 ]. 

In 1979 Lebeck [ 70,71 ] presented a mixed friction 
model including thermal rotation from frictional heating. 
Thermal and pressure caused rotation greatly alter face 
seal performance in the short run whereas after a period 
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of steady operation# the seal faces were parallel 
regardless of initial rotation. Thus this model can be 
applied to determine the effects of pressure and thermal 
rotation from frictional heating on performance of a 
face seal. In the same year Etsion and Etsion et.al 
[ 72-78 ] analysed the effects of coning misalignment 
and squeezing on the performance of a radial face seal 
considering both hydrodynamic and hydrostatic 
situation. . 

A strong coupling was found between angular 
misaligxment and transverse moment in both low and high 
pressure seals. Such coupling is a possible source of 
dynamic instability and can result in a wobbling of the 
primary seal ring. Hence# it can be concluded that the 
hydrodynamic transverse moment# which was overlooked in 
the past can affect mainly the dynamic stability of radial 
face seals and hence, has to be considered in any 
dynamic analysis of these seals. Further# the interaction 
between coning and diametral tilt has profound influence 
on seal performance. The axial separating force can be 
either positive or negative depending on the relation 
between the angle of tilt and the angle of coning. The 
same is true for tilting moment which can be either 
restoring or nonrestoring# and for the axial and angular 
stiffness of the sealing gap. Two different modes of 
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operation exist. In the first mode# contact between the 
mating faces# at maximum tilt, occurs on the outer diameter 
of the primary seal ring. In this case# high pressure on 
the inner diameter is needed for the stable operation. In 
the second mode, contact at maximum tilt is on the inner 
diameter of the seal ring, and high pressure on the outer 
diameter is needed for stable operation- Ruskell i 79 ] 
found that with suitable geometry good lubrication for 
face seal is possible. He carried out experiments also 
for describing the lubrication of these seals. 

In 1980 , Lebeck [ 80 ] and Young [ 81 ] brought 
out their experimental results showing how surface roughness 
and wear influence; load, seal-profile and leakage. Etsion 
et.al [ 82 ] analysed the effects of motion in 3 degrees 
of freedom of flexibly mounted seal ring on its stability 
and performance. A condition for angular stability is 
derived relating seal -operation-condition to its geometry 
and other design parameters. Walowit and Pinkus [ 83 ] 
gave a quantitative basis for the design and optimization 
of shrouded pocket face seals. It has been shown that 
for low seal pressures both hydrodynamic and squeeze 
film effects are important. 

Suyanami and Masuda [ 84 ], in 1982 , studied the 
behaviour of seal-ring in the shaft seal which operates 
on oil film to seal high pressure gas. The experiment was 
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performed on the 60 cm diameter seal ring. The analytical 
model is fojrmulated and predicts the oval deformation of 
seal ring due to hot spots and that the oval shape rotates 
slowly but periodically according to the shift of hot 
spots. Theoretical predictions are confiirmed to agree 
v/ith experimental results. Krauter [ 85 ] described an 
experimental programme concerned with the elastohydrodynamic 
behaviour of sliding elastomeric seals for the Stirling 
engine. During this programme an. experimental appjaratus 
was designed, built and used to measure oil film thickness 
distribution. For this oil film thickness measurement 
an optical interf erometeric procedure was developed. 

Etsion [ 86 ] analysed the Dynamic behaviour of a 
noncontacting face seal for the case of a rigidly mounted 
rotating seat and a flexibly mounted stationary ring 
taking into account various design parameters and operating 
conditions . The equations involved are non linear and 
solved mmerically. The effects of various parameters 
(for instance, radii ratio and coning) on seal stability 
is discussed and an expression for critical stability 
is formulated. 

1.5 PRESENT WORK 

At one time it was commonly believed that leakage 
flow across the faces of a mechanical seal was prevented 
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by intimate contact betv/een the two surfaces. Consequently 
seal design was oriented towards ensuring that face contact 
takes place over the whole face area and that wear is 
minimal under dry rubbing conditions which is achieved by 
r-sn dering the feces as flat as possible. In addition to 
this one face was flexibly mounted in order to follov; the 
inevitable run-out of the other. Majority of seals also 
have one of the faces made of carbon# which shows a consid- 
erable reluctance to seizure when in sliding contact 
with most engineering materials. Seal designs incorporating 
these features performed more or less satisfactorily 
for decades [39 ]. Hov/ever, in recent years duty 
requirements have grown more severe such as higher pressure 
and temperature leading to face distortions# 

It is now well established that for reliable 
operation over a long period, a radial-face mechanical 
seal must operate with its faces separated by a film 
of fluids [ 38 ], But the present day requirements 
necessitate face seals to function under different 
conditions . Seal leakage and short seal life are frequently 
observed in the equipments associated with seals. Sudden 
seal failure has made mechanical face seals most unpredict- 
able element. Thus in recent years a host of aspects 
influencing the operation of face seals have been studied. 
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Surface waviness, thermal distortions, surface roughness, 
vibrations etc. have been amply studied but we still fall 
short of adequate kno^;ledge that may provide a better 
seal design - 

Both hydrodynamic and hydrostatic effects were 
shown to have enormous influence over various seal 
characteristics [74,65 ] . Misalignments of the faces_ 
and coning are no less important in this regards [73 ]> 
Also, radial forces and squeeze effects are reported 
to be important in producing effects on seal performance 
[ 76,77 ]. In addition to this non~Newtonian effects 
do exist. Although non-Nev/tonian effects were reported 
to be insignificant [ 10,11 ], yet it is expected that 
a study of these effects when coupled with surface 
misalignment and flexibility of faces may lead to additional 
information that may help design a seal in a better way. 

The present work is designed to study some of 
these aspects of radial face seals. The seal faces 
are assumed to be narrow in that it renders Reynolds 
equations worth yielding an analytical solution in 
some cap*"*®, a consequence of elimination of circumfer- 
ential pressure gradient from the original Reynolds 
equation. In fact circumferential pressure gradient 
becomes insignificant when compared with radial pressure 
gradient [38 ]. Also this assumption renders the 
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curvature effects insignificant, as has been pointed 
out in [ 64 ]. Thus the present work may be divided 
into two parts ~ Newtonian effects and non-Newtonian 
effects. These effects have been coupled with porosity, 
flexibility and squeeze effects. Misalignment and coning 
have also been studied simultaneously. Centrifugal 
inertial effects and effects of radial forces have also 
been discussed. Most of the discussion is centered on 
the axial forces and tilting moments. These are really 
vital for the stability and durability of face seals. 

The Schematic of ®nd-face seal is given in 
figure 1.1 . 


f 



24 


CHrt-PTER II 

CENTRIFUGA.li INERTXA. EF FECT IN A &i.BAhljEli FACE SEAL 
2 . 1 INTRODUCTION 

It is now well established that, for a reliable 
operation of a radial face seal over a long period, it must 
be operated with its faces properly separated by a film 
of fluid. A proper gap is required to be maintained beca- 
use too small a gap will enhance the wear of the faces 
and too big a gap will produce massive leakage. In both 
the cases the problem of early seal failure is encountered. 
Also significant is the dynamic instability which is 
w:^ought by the vibrational motion of the primary seal 
ring. Dynamic instability has been observed experimentally by 
Haardt and Godet [ 50 ] and Etsion and Burton [75 ]. 
Investigators have considered several factors that may 
influence the performance of a face seal which continues 
to be an unpredictable element in the absaice of 
sufficient information. Lohou and Godet [ 41 3, Lebeck 
et.al [ 60 ] and Lebeck [ 61,62 ] and Etsion [ 72-78 ] 
considered the effects of misalignment. Surface waviness 
[ 42 3, thermal effect [ 54,55,57 3 » elastohydrodynamic 
effect [ 53,59 3# surface roughness [ 80,81 3 ^md squeeze 
effect [ 76 3 have already been studied, but there is 
no general theory as yet which can explain the strange 
behaviour Of a radial face seal. 
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Primary seal ring, being flexibly mounted on a 
support, remains sensitive to a very small variation in 
the operational conditions. The secondary seal, on accoxint 
of its high speed, may generate unavoidable amount of 
centrifugal inertia which in turn may well influence 
the performance of a face seal. 

Therefore in this chapter x-te start with a simple 
parallel face seal and analyse the variation of axial force, 
transverse moment and leakage after including the centrif- 
ugal inertia term in the equation of motion. 


2 ,2 MATHEMA.TICAL ANALYSIS 


Figure 2.1 is the configuration of a parallel 
face seal. It is considered that the flow is laminar 
and axisymmetric, the fluid is incompressible and the 
viscosity is constant. Further, in accordance with the 
theory of lubrication, all inertia terms have been neglected 
except that of centrifugal inertia. It is also assumed 
that the pressure is independent of the z-coordinate. 


The equations of motion, including the term of 
centrifugal inertia in it, are 


|E = p 


and 


1 


dz 


( 2 . 1 ) 


) 


(2 ,2 ) 
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where M is the viscosity P is the density and (u, v) the 
components of velocity i^ directions respectively, 
p is the pressure at (r/0). 


The equation of continuity in the cylindrical 
coordinates is 


I 1 - 

r ar 


(rv) + 


1 ^ 
r ae 


+ 


aw 

az 


0 


(2 .3 ) 


where w is the velocity in the z direction. 

The boundary conditions for the velocity components 
are ( o) is the angular velocity of the rotor). 

V = 0, u = 0 at z = 0 (2 .4) 

f V = O, u = rco at z = h (2.5) 

The solutions for u and v are obtained by solving 
equations (2*1) and (2.2), using boundary conditions in 
equations (2.4) and (2.5). Finally we get 


u = 


r ci)Z 

h 


( 2 . 6 ) 


and 

V z(z-h) - ^^tf-y (z^-zh^ ) (2 .7 ) 

"2M dr j2Uh 

Integrating the equation of continuity (2.3) across the 
film using the boundary conditions 

w =s O at z = O (2,0) 

w = 0atz=h (2.9 



we obtain 


i r V dO + i u d3) = (2 .10) 

z=h 

where L w represents the value of w at z = 0 and z = h. 

On further simplification, equation (2 .j9) reduces to 


1 

r 


a 

ar 


r { 


12^1 dr ^ 12 0 


2 t _3 
o hr 

U- 


1 ] = O 


(2 . 11 ) 


Equation (2.lJ|)) is the final equation governing the flow 
of the fluid film between the faces of seal. 

Since h = constant, equation (2.lfi) reduces to 

a_[ Hg - Pu?r). ] = o C 


Equation (2.1^) is solved with the help of the following 
boundary conditions 


p=p^ at r=rj^ and p=P 2 at r = r 2 


(2 . 13 ) 


where r^^ and r 2 are inner and outer radii respectively. 
Pj^ is the inner pressure and is the outer pressure. 
Thus the pressure distribution is given by the following 
expression 


Llog(| — )(p^-P2 )( — 

1 logcg.) 


-{|5-Puh4-r2)i{ log(|^) 


log(~) 

^2 


+ Pi ) J] 


(2 , 14 ) 
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2 *2 .1 Axial Force I 


The axial force is defined by 


p 27T r dr 


( 2 . 15 ) 


On substituting for p form equation (2.14) in equation (2.15) 
and integrating, we get the exoression for F as boloxv 


= I [£ 2^1 log (^)-(r2-r2)HCpj^-P2) 

3 A 2, 2 _2s, , 1 

- 20- Pu ( 


log(-i) 


+2pj^(rj-r^) - Pa^CrJ-rj) ] (2.16) 


2.2.2 Restoring Moment t The tilting moment about x axis, 

V 

also called restoring moment, is defined by 


^2 2 

M =-/ /^pr cos© d© dr 
o r^ 


(2 .17) 


2 2 

The factor pr in the integrand pr cos© in eeguation 
(2.17) is a function of r only. Further, when cos© is 
integrated between the limits 0 and 2n, the integral becomes 
zero. Therefore 


M = 0 
X 


(2.18) 


2.2.3 Transverse Moment : 


The moment about y axis is called the transverse 
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moment and is defined by 

$TC ^ 

M = / / 2 pr" sine de dr (2.19) 

^ or 

Substituting for p from equation (2.14) in equation (2.19) 
and integrating, the expression for becomes 




Y 


JT 

2 [~ I { log (^)}{ 


log(-~) 

2 


1 / 3 3 \ r / \ 3 2 / 2 2 \ ^ 


+ 3Pi(^2'‘^l^ ’lOO 3r2+2rj^-5r^r2) ] (2.20) 


2*2.4 Leakage : 


The leakage is the flow across the seal ring and is 
defined as 

2 Hr 

Q = S J V r dQ (2.21) 

o ^ 

Substituting for v from equation (2.7) in equation (2.21) 
and integrating we get 


Q 


^ iL 

6 M 


3 


[{ 


log(— T") 
2 


ifCpi-Pj) 


3 

120 


2,2 2 ^ , 

Pg) (r,-r„)} 


( 2 . 22 ) 


2.3 NON-DIMENSIONAL SCHEME 


We adopt the following non-dimensional scheme 
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"2 

2F 


' Pi = P2 ' 


A = 


r 2 2 

Si 


' ^1 ~ ' 


(2 .2 3 ) 


F = 

z 


^P2^2 


M = 

y 


; Q = 


Q 


^2^2 


( Tth /6p^M) 
The non-dimensional axial force becomes 


F. 


L {2 logC—r*)-^ l-?3_ — ;;; — }{(pj_-l)~ 


3 


•1 


log(r^ ) 


^ -if ] (2 .2 4) 

The non-dimensional form of the transverse moment is I 


+ 2p^(l-rj, ) - 20"'' -1 


M = [{ jlogC^ — ) -|■(l-r^)}^ — }x { (p^-1) 


log( r, ) 


~ 20 ^ 3 Pl^^“^l^ 


1 ;; n =3 
3 Pl^^' 

1 -Y / ^ r5 _ —2 


100 


A (3 + 2r - 5r^) ] 


(2 .2 5 ) 


The dimensionless form of the side leakage is I 


Q = - { 


log( r, ) 


} L ( P^-l ) - lo (2.2 6) 


2.4 RESULTS AND DISCUSSION 

The dimensionless axial force the transverse 

moment Hy and the leakage Q are calculated from: equations 
(2 .2 4)/ (2 .2 5) and (2.2 6) respectively. The results are 
presented in the form of graphs in the figures 2.2/ 2.3 and 


2 .4 
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iiach of the quantities F , M and Q is a function 

z y 

n 2 2 

Po r- 

of the inertia parameter a = and the radii ratio 

r^ = p— . Kf for a fixed value of r 2 # and the density P/ 
would depend on the angular velocity (*> . Thus the minimum 
value of X can be zero v/hich corresponds to the hydrostatic 
situation. The maximum value of X would depend on the 
value of 0 ). In the present chapter we have taken the 


maximum value of X to be 350. It may be even greater. 


Further , the radii ratio r^^ would increase or decrease 
according as r^ increases or decreases, or r 2 decreases 
or increases. In each case increasing r^ would signify 
decreasing the thickness of the seal ring. The present 
analysis assumes the value of r^ between .80 and .95 v/hich 
is in accordance with the real situation. 


Figure 2 .2 represents the variation of dimensionless 
axial force F^ with the inertia parameter A. This figure 
shows that, for all values of radii ratio r^^, the dimen- 
sionless axial force F increases as X increases* However, 

F decreases with an increase in the value of r, . The . 
z i 

increase in the value of F„ with X is sharper for smaller 

\ F / Z 

value of rj^ than for large value of r^^. 

Figure 2 .3 represents the graph between the non- 
dimensional transverse moment for various values of 
r^. This figure shov/s that the transverse moment 
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decreases with an increase in the value of both A. and r^. 

However, the rate of decline in the values of M vrith 

Y 

an increase in A is less for larger values of r^^ than 
for smaller values of r^^. It tends to become almost 
constant with respect to A as r^ approaches' ‘unity* 

Figure 2,4 represents the variation of non-dimension- 
al leakage Q with A and indicates that Q increases with 
an increase in both A and radii ratio r^. 

The aforesaid discussion leads to the conclusion 
that for a fixed value of r^, higher angular velocity 
results in higher values of axial force and leakage* 

Thus unless additional devices are used, the lubrication 
system in which the rotor has very low angular velocity, 
would witness an early failure of face seal because 
of massive wear in the abs ence of adequate separating 
force. On the other hand the systans in which rotor 
has very high angular velocity will starve for lubricant 
because of the massive leakage it will have to undergo. 

Thus a proper angular velocity of rotor is required to be 
maintained to ^^^^re smooth f-unctioning of face seal 
for a long period. 






Fiq 2 3 Dimensionless transverse moment as a function of inertia 
parameter for various values of radii ratio. . 




i; ■ ir,iriid‘ p^romct*rS 
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CHAPTER III 

SQUEEZE EFFECTS IN J^IISALIGNED RADIAL FACE SEALS 
WITH CONING 

3 . 1 INTRODUCTION 

The primary seal ring of an end-face seal is 
flexibly mounted on spring or bellows (Fig.l.l). 
Imperfections in the axisymmetry of such a support 
prevents the primary seal ring from being ideally 
perpendicular to the shaft. Since the clearance between 
primary real ring and seal seat is of the order of few 
micrometers, and since imperfections of this order of 
magnitude are almost inevitable in the flexible support^ 
it is very likely that diametral tilt is inherent in 
end-face seals. 

Coning of the mating faces is also inevitable 
due to thermal and mechanical distortions . Hence to 
predict accurately end-face seals' performance, both 
diametral tilt and coning have to be considered. 

There is considerable scatter in the life of 
radial face seals. In an attempt to visuatize the crux 
of the problem, Sharoni and Etsion L 72 ] considered 
the effects of coning on seal characteristics, Etsion [76] 

This chapter has been published in the form of following 
research paper "Squeeze Effect in Misaligned Radial Face 
Seal With Coning", Wear, Vol,85#1983,p-l43« 
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analysed radial face seals by treating the film 
thickness as time dependent and considering only misalign 
ment of seal faces . 

The analysis of ref. 72 has indicated that high 
pressure seal stability may be improved by coning. In 
the present chapter the squeeze effects in a misaligned 
high pressure radial face seal with coning is analysed. 
The film thickness was regarded as time-dependent and 
expressions for forces and moments have been obtained. 
Etsion's [ 76 ] analysis may be obtained as a special 
case of the present analysis. 


3.2 MATHEMATICAL ANALYSIS 

The Reynolds equation for a misaligned radial 
face seal (fig. 3.1) using an incompressible fluid is 



(rh^ -1^) 


ar 


^ 1 -^(h^ 

r ae ae 


6Mr(m + 2 

00 ot 


where the film thickness h is (figure 3.1) 


(3.1) 


h = c + r r cos© + P(r-rj^) (3.2 ) 

where c is the seal clearance along the central line 
joining the centres of seal faces, r is the angle of tilt 
and (r,©) is the coordinate of the point through which 
the thickness is measured, P is the angle of coning and 
rj^ is the inner radius of the seal ring. 



40 


Equation (3.1) may be solved separately for pure 
hydrodynamic/ pure hydrostatic or pure squeeze effects 
but only the squeeze effects are considered here. Using 
the narrow seal approximation [ 64 ]/ equation (3.1) 
reduces to 


i__(h3 ^) _ 12M ^ 
ar ^ ar at 


(3 .3 ) 


where 


ah 

at 


is obtained by differentiating the film thickness 


h in the equation (3.2 ). Thus 


Vo 


at 


-v^+ 7 r cos 6 

c m 


(3.4) 


where is the velocity at the centre, i.e. v 


at 


and f 


AJL 

at 


The symbol r^^ denotes the mid radius 


i.e. r = ( — o— ^) where r^ is the outer radius of the 
in ^ ^ 

seal ring. By combining equations (3.3) and (3.4) we 
obtain 


-^) = 12M(-v+t r cose ) C3.5) 

ar ar c 

To integrate equation (3.5) the pressure boundary conditions 
are 

p = O at r = rj^ and at r = r 2 

-|§- = 0 at some r = r' (3.6) 


where r' is the local radius corresponding to the 
pressure extremum. 
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Integration of equation (3.5) with respect to r gives 


h^ ■§§“ ~ 12M(-v^t f r^ coseJr + K^C©) (3.7) 

where K^(0) is the constant of integration which is a 
function of 0 . To evaluate 0 we make use of the second 
pressure condition given in the equation (3.6). Thus the 
value of K^(9) is obtained as 

K, (e) = -12M(-v + r r cose )r' (3.8) 

i 0^1 

On combining equation (3.7) and equation (3.8) we obtain 
final expression for the pressure gradient as the following 


^ = 12M 


/ “V _+ f r cosS ) ( r-r' ) 
\ ^ rn 



(3.9) 


Integrating equation (3.9) we get 

p(r) = 12M(“V + ‘f r cos 9^ -- - i- — ^ dr +KL(0) (3.10) 

h^ ^ 

where ^^(6) is the constant of integration which is a 
function of 9 , Equation (3.10) and equation (3.6) 
together give K 2 (9 ) =0 and 



r dr 




(3.11) 


By solving equation (3.11) and then substituting the value 
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of r* in ecmation ( 3 .IO) and carrying out some algebraic 
manipulation we get the expression for the pressure as 


ehmCv-i r cos0 ) (r-r, ) (r^ -r } 

c in ± £ 


{ c+7 r COS0 +P(r-r.,)} {c+'X r^ cosG+P (r^-r^ )} 


"ra "1' 


( 3 .12 ) 


3.2 2. Axial Force 


The axial force is defined as [ 76 ] 




^2=2 r^f S ^ P dr de 
■°^1 

3 .2 .2 Restoring Moment I 


(3.13 ) 


The restoring moment about x axis is defined as 


M 


2 t To 

~2 r^ J* J ^ COS0 p dr d0 
0 

3 . 2.3 Translational Damping Coefficient I 

is a function of both v^and f as does 
The translational damping coefficient is defined 
as the value of axial force F_ when we put t = 0 in the 
expression of F_. Thus 


( 3 . 14 ) 


P = F^ (f = O) 
V z 


( 3 .15 ) 


3 . 2.4 Cross -coupled Damping Coefficient I 

The cross-coupled damping coefficient F^ is 
defined as the value of axial force F when we put v = O 


in the e^qpression of F . Thus 




F^(v = O) 
z C 


(3.16) 
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3.2.5 Rotational Damping Coefficient: I 

The rotational damping coefficient is defined 

as the value of the restoring moment M when we put v = 0 

X c 

in the expression for M - Thus 

M. = M (v = 0) (3.17) 

7 X C 

3.3 NON-DIMENSIONAL SCHEME 


The following sxibstitutions are used 


r . - ^1 . - ^ 

^ ^ ^ = HT ' ^m = r. 


rr- 
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V 
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Using this scheme to obtain the dimensionless pressure 
as the following 


P 


6 M( — )^ ( — - 7 ' it cos6 ) (r-r. )(l-r) 

C XV) HI X 

{ 1 +er COS 0 + 6(r-r2_)i £L4€ r^^ cose + 


(3.19) 


Now let us suppose 

1 ^ 1+ecose +6(l-r-i ) 

I = / p dr = rr [log { } 

« (ecosG) 1+e r, cose 

(l-r, ) (ecos 6 +6 } 
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m ml 


(3«20) 



Using equation (3.2 0) ‘in the formula for the translational 

} 

^damping coefficient the rotational damping coefficient 
and the cross-coupled damping coefficient , we 
obtain the non-dimensional integrals for their expressions 
as the following 

V m *' 

o 

_ m ■! r 1+ecose +6(l-r, ) 

= '3 [log £ i-i 

0 (ecos0) l+er^cosS 

(1-r, ) (ecos 0 +5 ) 

-{ i } ] de (3.21) 

1+er cos8 +6(r ) 

m ml 

and 


F . = M 
r V 


- 2 1 
= -r / 
m *' 


o 3 2. 

e cos 6 


[log { 


1+ecos 9 +6 (l-r, ) 


1+er^cose 


-£ 


(l-r^)(ecose +6) 
1+ ^r^^cose 


■} ] de 


( 3.22 ) 


and 


- t 2 

= - rl / I cos 0 de 
7 m«' 

o 


-3 

= 


n ^3 
O £ cos 


[log{ 


1+ecos 0 +6(l-r, ) 


■} 


l+er^cose 


(l-rj^)(ecose +6) 


} 3 de 


(3.23 ) 
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3.4 RESULTS AND DISCUSSION 

The dimensionless translational damping coefficient 
the cross-coupled damping coefficient F^ or and 
the dimensionless rotational damping coefficient are 
evaluated by numerical integration of equations (3.2 1), 
(3.22) and (3.23) respectively. Each of F^/ F^ is a 

function of tilt parameter e and coning parameter 6 . The 

minimum value for both £ = ■ " and 6 = is equal to zero. 

This corresponds to the cases of no angular misalignment 
and no coning respectively. The expression for the film 
thickness from equation (3 .2 ) can be written as 


— = l+£rcos0 +5(r-r, ) 
c 1 

Q = % corresponds to minimum film thickness (figure 3.1). 

Thus when 0 = n: , the film thickness, because of the 

misalignment, is minimiom only if e 1* But if £ > 1 #_a 

negative film thickness is obtained which is meaningless; 

Thus the maximiom value that e can take is 1. Further, 

when r = r 2 / 6 = and e = 1, the non-dimensional film 

thickness - = 6(l-r, ). For all practical purposes the 
o 

value of r^ can be as near to 1 as .98 (this means the 
thickness of seal ring is very small). Thus when = .98, 

= 6 X .02 . Further, ~ < 1 at the end where 9 
and r = r 2 . Therefore, the maximum value that 6 can assume 





is the one which keeps the value of ~ = 6(l-r^) smaller than 
1. Which is possible when 6 < 10* Thus in present 
discussion the maximum value of 6 is taken to be 10. 

Each of the characteristics has been plotted as 
a function of the tilt parameter t, the value of which 
vary between 0 and 1, for the values of 6 varying betv/een 
O and 10. 

The positive values of the coefficients and 
indicate that if the seal ring is disturbed from a 
certain position (c,t) , the response and restrict 
the disturbance. 

The cross-coupled damping coefficient shows a tendency 
to restore any given tilt of the faces, that is a positive 
f which increases the tilt angle O', generates a separating 
force F^ this in turn increases the clearance c. 1- positive 
V which decreases the clearance c, produces a moment 
which reduces the tilt O'. In both cases, the result is 
a stabilizing effect tending to prevent any contact between 
mating faces. 

Figures 3.2 and 3.3 indicate that the translational 
damping coefficient remains almost constant at small tilts 
and increases sharply as the tilt parameter approaches 
linity. 



Figures 3.4 and 3.5 show that the cross-coupled 
damping coefficient depends strongly on the tilt parameter. 
At small values of tilt parameter e, which results from 
Small tilt angle Y, or large seal clearance C/ the coupling 
effects are small. As e increases and approaches iinity/ 
the cross-couple damping coefficient registers a sharp 
increase. 

Figure 3.6 and 3>7 indicate that the behaviour of 
the rotational damping coefficient is essentially the 
same as that of the translational damping coefficient. 

All of the graphs show that as 6 increases, the 
values of forces and moments decrease. They decrease 
with an increase in the value of radius ratio. 

Finally, it may be concluded that squeeze with 
coning may be desirable where the large values of the 
separating forces produce leakage, because coning 
causes a reduction in these forces. 






^ 0-0 0.2 0.4 0-6 0.8 10 

6 — » 

Fig 3 5 Dimensionless cross -coupled damping coefficient as a 
function of tilt parameter for various volues of coning 
pai ameter 




Fig. 3.6 Dimensionless rotational damping coefficient as a fi 
of tilt parameter for various values of coning para 
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CHAPTER IV 

HYDRODYISiAMIC EFFECTS IH POROUS RADIAL FACE SEALS V7ITH 

MISALIGNMENT 


4.1 INTRODUCTION 

A radial face seal (fig.4*l) is probably the most 
non-predictable machine element. However, face seals 
are widely used in a host of industrial equipment and 
in other applications. The inability to predict seal 
life within reasonable accuracy is a major disadvantage 
and seal failure is a common problem. Recently many 
investigators have treated the problem and there have 
been several hypotheses put forth to explain the mechanism 
responsible for the development of lubricating film pressure 
that acts to separate the primary seal faces. These 
hypotheses include surface angular misalignment L 28,29], 
surface waviness [ 2 7,42 ], microasperities [ 30 ] and 
thermal deformation [ 57,58 ]. Inspite of these investigations 
the theory of face seal lubrication demands still more 
attention. 

The restoring moment which is an imp03rtant factor 
for stabilising or destabilising the seal-function has 
been discussed in [ 52,55,65 ]. The transverse moment 
that leads the angular misalignment vector by 90*^ is 
generated by hydrod 3 mamic effects [ 74 ]. Tlie author 




56 


has shown tha-t the transverse moment due to the angular 
misalignment of the seal faces can result in a wobbling 
of the primary seal ring. It is fiirther concluded that 
the hydrodynamic transverse moment certainly plays an 
important role in radial face seal operation# 

The analysis made by [74 ] include only angular 
misalignment. However, a face seal may be sufficiently 
influenced by several other factors depending on the 
operational conditions and nature of the face material 
of a face seal. This is significant because, apart 
from influencing other seal characteristics, these 
may influence transverse moment also which determines 
dynamic stability of a face seal. One of several factors 
is the porosity of the seal faces. 

The porosity of the faces, which has not yet 
been analysed for misaligned face seals, constitute the 
main objective of this chapter. The primary seal ring 
is considered to be coated with a porous material 
(figure 4.1) and the fact-s are considered to be 
misaligned. 

4.2 MATHEMATICAL ANALYSIS 

The Reynolds equation for a narrow seal with 
incompressible fluid as lubricant is * 
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'.'.■'I'lCiiTG h is the film thiclcness/ Vj^ is the velocity of 
fluid at the face of secondary seal at 2 =r h and is 
the velocity of the fluid at 2 =0 i,e. at the face 
of primary seal ring which is porous, p is the pressure 
of the fluid between the faces at a pt. (r,0 ). 

Assuming that the flow in the porous shell in 
fig. 4.1 is governed by Darcy's law L 87 ] 




2 


1 3P 
r 5B 




2 ^ 
M dr 


and 

Q( 3P 

"Jz = ” M 3z~ 

where ^ is the permeability, M is the viscosity of 
the fluid P is the pressure in the porous medium 
and q^, are the respective flux of the fluid in 
the 9 ,r and z directions. 


(4.2) 
(4.3 ) 


(4.4) 


The equation of continuity in cylindrical coordinates 
system is * 




39 




Equation of continuity may also be written as « 
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+ 


rae 


*<qn ) + 


az 




o 


(4.5 } 


Since the circumferential pressure gradient in 
annular cases are much smaller than other pressure 
gradients [38 ], the terra involving can be 
neglected in the equation of continuity ( 4 . 5 ) 
because i- is neglegible as compared to • 

Thus on putting the values of q^ and q^ from equations 
( 4 . 3 ) and ( 4 . 4 ) in equation (4.5), the equation of 
continuity reduces to 


3 ^ 3 ^ 

ar^ dz^ ~ ° 

({?■ and M are independent of r and z). 


(4.6 ) 


dP 

It is reasonable to assume that gg— is linear across 
the porous matrix [87 }, Finally from equation (4.4) 

ap 


we have = 0 at the outer surface of the porous 

shell, since there is no axial flow of the lubricant 
there. 


These assumptions can be expressed as 

= ^(z+H) “and ~ — I H = O ( 4 . 7 ) 

az 3z lz=:-H 

where ^ is a constant depending on r only and H is 
the thiclcness of the porous shell, i.e., z = -H refers 
to a point on the outer surface and 



(4.3 ) 



59 


From equation (4.7) we have 


3%> 


az 




1 

(z+H) dz 


(4.9) 


On using equation of continuity (4.6) here, we obtain 


az ~ ^2+H) — ^ 


a% 


= -(z+H) 

3z^ ar 


(4.10) 


Furthermore^ we have =0 at z r= h because secondary 
real seat is not supposed to be porous# Also# the 
pressure is continuous at the interface ^i.e# prrPatz = 0. 

In view of actuations (4.8) and (4.10) we have 


Q( 

V r: ^ H 
O /i 




(4.11) 


Also, Equations (4.6) and (4.9) together give 


a^ 

17 


a^ 


az 


_ ^ 


( 4.12 ) 


On using equation (4.11) in equation (4.1) we have the 
final governing equation as I 

^) = 6^^o)r 120Hr -^-1- 

ar ar 30 0r* 


(4.13 ) 


Using the analysis for narrow seals in [ 64 ], equation 
(4 . 13 ) further reduces to the following 

|^LCh^+12pH) ||~] = 6W0 


ar 


ae 


(4.14) 



60 


V7here h, the film thickness, is expressed as 
h = c + T r COS0 

where c is the clearance of the seal faces along the 
line joining their centres and r is the angle of tilt. 


Now 


= -r r sine 

Substituting for from equation (4.15) in equation 

(4.14), we obtain 

“gj?' L (h^+12 0H) = -6M 0 ) yr sin0 

Integrating equation (4.16) with respect to r we get 


(4.15) 


dp 

~W 


-3M 0 ) rr sine + K(e) 

(h^tl2(?H) 


(4.16) 


(4.17) 


where ic(e) is the constant of integration which is a 
function of © . 

Therefore 

r [(-3M u yr‘ sine ) + K(e) ] 

p(j-) _ f <ar 

(h'^+i2orH) 

We have the following boundary conditions I 
p(r2) = 0 = p(r2 ) 

Therefore 

To (~3M 0 ) 7x2 + K(e) ^ 

p(r2 ) = O = /^ [ -5 ; 3 <ar 

2 r (h^+129H) 


(4.18) 
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which yields the value of K(e ) as I 

r 2 

/ ^ 3 A 0 ) Y r sin ft 

(h^+ 12 pH) 

K<B ) = — (4-19) 

f^2 ^ 

(h^+12pM) 

In view of [ 74 ] , in the absence of hydrostatic 

component, cavitation is assumed to occur in the region 

^ ah ^ 

where 0 i.e. v/hen ^>71 and thus pressure is asstmed 

to be zero for 6 > 7 i . But when the hydrostatic component 

does exist, the extent of cavitation is dependent on the 
pressure differential across the seal boundaries. The 
cavitation zone decreases as the sealed pressxire increases 
until a full fluid film condition is reached. We shall 
now consider the cavitation flow. 


4.3 CAVITATION FLOW 
4.3.1 Axial Force I 

In the case of very small pressure differential 
across the seal boundaries, the axial force is defined 
as [ 74 ] . 


z 



pr dr d 0 


Integrating the right hand side of the equation by 
parts we get a more simple form for F^ as the following'. 
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^2 ^ ^ / ^ ^ ^ ( 4-2 0 ) 

On putting the value of — from equation (4.17) in 
equation (4.2 0) 

1 r'^^^2! “3^ o 7 sine + K(e ) ^ 

~ TJ f L ■ IT ] r dr de ( 4 . 21 ) 

o (c+r r cose ) +12pH) 

4.3.2 Transverse Moment I - 


The moment about y-axis is given by 


^ ^2 2 

M = / / p r sine dr de 
7 o 

Integrating this with respect to r, by parts, and then using 
equation (4.17 ) we get 


, Ti r„ ^ -3M o r r sine + K(0) _ 

= - i / / 2 [ 3 ] 

^ ° ^ cose)"’ + 12 pw) 


sine dr d0 


4.3.3 Restoring Moment I- 

The restoring moment is the moment about x-axis 
and is given by 

r^2 2 

r ^ pi 

o r- 


#v ^ O ^ 

M = - / J ^ p r cose dr d6 


X 


■1 

= T //^ [ - 3 ^ “ - 1-^(3 ) ]^2 

(c+7 r COS0 )"^ +12 011 
4.4. NON-DIMENSIONAL SCHEME 


(4.22 ) 


(4.23 ) 


Various non-dimensional quantities are defined 


below 
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r r- 


£ = 




a 


; p = 2=. 


M w r: 


12 OH . - 

, F 


M 6. J ' ’"y 


; K(e) = mi^ 

M £i) 

2 

M c 

-JZ 


M £i) r;: 


M 


M C 
X 


X 


25_. 


-2 


M 6) 


s ' ar 


M o 


02 _ -3£ r sin ft -»• K(9) 

(1 + e r GOS0 )^ + a 

(4.2 4) 


Thus K(0 ) = 


/ 


J 


3£ sin 9 r dr 


[d+e r COS0 + a ] 


or 


•1 (1+e r cos 9)^ + a 


(4.25 ) 


and 


5(f) = / [ (-3e f sine^. K(e) ] ^ 


■1 


(1+e r COS0 + a 


The non-dimensional axial force becomes 


F = 
z 


F c 
z 


(4.26) 


M Cl) r^ 


= i// [ as . / 3 ,. ine _ ^ (e2 ]j2 


p J J ^ *1 — 

o (1+e r cos9)'^+ a 


de 


(4.27 ) 


The non-dimensional fom of the transverse moment 
is given as below 


1 f /I [ 3£ sine-%.i,3 j3 


M 


^ 3 - (1+e r cos0)^ + a 


( 4.2 8 ) 
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The non-dimensional form of the restoring moment 
becomes 

M = i / r-3£ r^ sine, ^.„K(.e) n -3 ^ (4.29) 

^ ^ o (l+e r co^ ) + cc 

Equations (4.27) ^ (4.2 8) and (4.2 9) when evaluated 

give the values of the seal charactaristics to be consioered 
for cavitation case. VJhen full fluid film prevails, 
the component of restoring moment becomes zero whereas 
that of the transverse moment is dovibled. 


4,5 REDUCTION OF DOUBLE INTEGRALS TO SINGLE INTEGRALS 

To find out an expression for p(r).’ let us asstime 

j f ^ t k( 9 ) (4.30) 

” _ (l+e r cose)^ + a 

Let = a, 1 + e r cose * u dr = g^cose ^ £ co^ 

(4.31) 


With these substitutions our integral I reduces to 

2 


3 sine [ r 


U 


2 ^3- 

e cos e 


+/ ^ au] (4.32 ) 

^3 + 0.3 + <x3 


where 


K(Q )£cos% 
3sin9 


On integrating (4.32 ) we obtain 

I ^ - 3slne [ + MI 2 ] 

0^2 e cos 0 
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where 


I, = a' 


log(u3+a^)+%2«i log ),il±4gj 

u^-au+a^ fa Vs a 


(4.34) 


and 


= log 


+.4- tan“^ (2 u-a ) 


( 4.35 ) 


Thus using equations (4.2 6 ), (4.33 ), (4.34 ), (4-35 ), we can 
evaluate p(r) easily for various values of a and e« 

Proceeding in the same way, we can reduce the double 
integrals in equation (4.2 5 ), (4.2 6 ) and (4.27) to single 
integrals by integrating these integrals w.r.t. r. Thus 
the single integral expressions for separating force, 
transversse moment and restoring moments are as the following 

(6-M) 


p _ ■ J tanS sec^d [ ^ 


26*^ o 


- 4u + 


log(u^+ot^ ) 


Qo 2 2 

+ Qj^log(u4a) + 2 — log(t -au+cc ) 


(2Q^+Q^cc) 

■^3 a 


tan 


-1 (2u-g ) 

Y3 a 


u- 


3 de 


(4.36) 




where Q-, 




b,oc-L- 

♦ r\ X n 


3 a 


3 a 


_ _(4«2M+a'^ ),*L^=(4a +1-M); u=l+e r cos8 ,* 


uj^ = 1 + e % cose ; uj = 1 + e cosG 



66 


And 


M = - 
X 


^ tan^e sec^e [ 7 ^ + 


z o 


( IO 4 M )u 


(lO+SM+a^ .7 7 

■ 3 log(u +a-^ ) + log(u -au-i-a ; 


H" Q log C ^ 

^ 1 ( 2 u-a) 


■f3 a "fSa 


tan 


T3 a 


] 


u. 


( 4 . 37 ) 


And 
M = 


1 j- ^ ~ 4;. 

5 ~J' tan0 sec 0 [ j— - + (10+M)u 

e o 

- U og(uV^ ) + ^log(u2-au+a2 ^ 

IfS a 


+Q^j^log(u+a )+ — ^ 


a 


-tan 


where Q 


11 


“ ' %2 - — 03 — 


(4.38) 


3 a 


3 a 


33 


^11°^-^2L22 . 

3CC ' 


= (5+3M+5a^) # L 22 = “ ^ • 


3.5 RESULTS AND DISCUSSION 

These are two expressions for each of the quantities 
separating force# restoring moment and transverse moment. 

In each case one expression involves double integration 
and the other involves single integration. Each of the 
expressions (equations (4.36) , (4.37) and (4.38)) involving 
single integration also carries a factor of some power 
of cos0 in the denominator of the integrand. Since the 
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integrand has to be integrated between the limits 0 = 0 and 
0 = Ji , the integrand becomes indefinitely large in the neigh- 
bourhood of and at the value of 0 = 7T/2 . Thus numerical 
integration of these expressions results in frequent m.assages 
of overflow. Thus double integration of equations (4-2 7 ), 

(4.2 8) and (4.2 9) is performed to obtain the n^americal values 
of the axial force and moments. 

Before performing numerical Integrations, a value of 

the porosity parameter a = — ^ is required. Various porous 

c 2 

materials have permeabilities ranging from 100 x 10 cm , 

to abqut 2000 x 10“^ cm^ - Thus for a bearing material 

having"^ 9 = 1500 x lO”^ cm^ , H = 0.508 cm and c = 2 .54 ic lO'^cm, 

is about 0.54. Thus the values of a = are taken 

3 c 

c 

between 0.2 5 and 1.0. 

The axial force F is presented in fig^lres 4.2 and 
4.3, the restoring moment in figures 4.4 and 4.5 and the 
transverse moment in figure 4.6 and 4.7. Each of tnese 
quantities is plotted against the values of e ranging between 
O and 1, for wide range of values of porosity parameter a. 

An increase in the value of « may either be caused by an 
increase in the value of permeability pr or by a decrease 
in the value of seal clearance- 

Figure 4.2 and 4.3 indicate that the value of 

axial force increases sharply for an increase in the 

.^111 it attains some value about 


value of tilt parameter e 
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.2 5. For e > .2 5, an increase in the value of e is 
followed by an increase in the value of axial force but 
less sharply. Figure 4.4 and 4.5 show that the value 
of restoring moment increases very sharply with e but 
the sharpness slackens as e approaches \inity. Transverse 
moments [ figs. 4.6 and 4.7 ] behaves essentially as the 
restoring moment does. 

Each of the figures indicate that each of these 
three quantities decrease with an increase in the value 
of a. Thus when there is a case of 'excessive leakage 
and it is sought to lessen the transverse moment to 
bring about stability in a face seal, the porosity of 
face material may serve the purpose, since it tends to 
decrease them. Also, when compared with ot =0, i.e., 
a case when no porosity is there, the porous face seal 
has a tendency to undergo less shooting incrments in the 
values of these characteristics with e . 

It has been established [ 74 ] that the hydrodynamic 
transverse moment which can affect the dynamic stability 
of radial face seals has to be considered in any dynamic 
analysis of seals. It is , therefore, concluded on 
the basis of present analysis that the porosity of face 
material has important influence on the stability of 
radial face seals because porosity as such significantly 
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influences the dynamic transverse moment. Further 
reducing the hydrodynamic transverse moment to a minimim 
may result in more stable seals. 








Fig 4 4 Dimensionless restoring moment os o function of tilt parameter 
tor vonous values of porosity parameter. 





o 
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' CHft,PTJ» V 

EIASTOHYDROSmTIC RADIAL AtqPi AXIAI. FORC3BS IN A MISALIGNED 

RADIAL FACE SEAL 


5.1 INTRODUCTION 

The secondary seal in a radial face seal is usually 
an elastic member such as an O ring or a piston ring. Its 
main task is to prevent leakage via the radial clearance 
between the primary seal ring and the seal housing. Since 
the primary seal ring is flexibly mounted on a support# 
it has five degrees of freedom L 86 ] out of which two 
are the two perpendicular radial directions* Further# 
since the leakage highly depends on the eccentricity between 
the seal seat and the primary seal ring [, 29#3l 3 # ^ shifting 
of primary seal ring in the radial direction is followed 
by a heavy leakage which in turn brings about seal failure. 

Both hydrodynamic and hydrostatic effects can 
produce radial forces in a misaligned radial face seal. 

The hydrodynamic ccxnponent is due to nonaxisyiwmetric 
tangential and radial shear and the hydrostatic component 
is mainly due to nonaxisymmetric hydrostatic pressure 
distribution L 65 ]• 

In order to provide a better -understanding of 
radial face seal mechanism Etsion and Sharoni [ 77 ] and 
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Etsion [ 78 ] have studied the effects of radial forces 
combined with the effects of coning on a radial face 
seal. It has been shown that radial forces have direct 
influence on the stability of the seal. Axial forces, 
which have impact on leakage and wear# have also been 
studied for various cases [ 64# 65»72 ,73#74#75#76# 86 ]. 

But these analyses ignored the consideration of the 
elastic nature of the seal faces which may be considerably 
deformed elastically xinder high pressure. As a mater 
of fact#the film thickness undergoes a change and 
accordingly the radial and axial forces follow suit. 

It is thus the objective of this chapter to study 
the effect of elasticity of seal faces by considering 
film thickness as a linear function of pressure. The 
influence on the radial force, which is so much vital 
for the stability of a seal, has been discussed in 
detail. In addition to this, the changes in the axial 
force have also been discussed. 

5 .2 MATHEMATICAL ANALYSIS 

Figure 4,1 depicts the geometry of a misaligned 
radial face seal. The Reynolds equation for a narrow 
seal with an incompressible lubricant is . 




ay r 


ae 


(5.1) 
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where is the film thickness when the pressure is p* 

In fact 

h' = h + sp (5.2 ) 

where h is the film thickness when no elastic effect 
is produced, h for a narrow seal is given by 

h = c + T r COS0 (5.3 ) 

m 

s is the parameter which is inversely proportional to 

the modulus of elasticity of the face material. Further, 

c is the clearance along the centres of the faces of a 

seal and Y and r are the angle of tilt and the mid 

m ^ 

radius respectively. 

For a solution which deals with the hydrostatic 
component only and for a narrow seal [ 64 ] equation (5.1) 
reduces to 


t- = ° 

The pressure boundary conditions are 


(5.4 ) 


p = p^ when r = r^^ 
p = when r = r 2 


(5-5 ) 


=s O at some r = r^ . 

Integrating equation (5-4) and using equation (5.2), we 


get 


yi (h+sp)^ ] = B 


(5.6) 
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Integrating this equation we get 

(h + sp)^ = Br + C (5.7 

where B and C are constants of integration depending 
only on ©. These constants are determined with the 
help of the boundary conditions in the equation (5.5). 

Thus we have 


B 


L (h+spj^ )'^-(h-t-sp^ 3"^] 


(ri-r^ 


} 


(5.8) 


and 


C 


Lr^(h+sp^ 


(ri-rj ) 


(5.9) 


Now for a narrow seal [ 64 equation (5.7) may be 
approximated as 


(h+sp ) = 


Br„+C 

m 


[ (h+sp^ )^-(htsr:^ )^]r^+[r^(h+sj:^ ^^”^2 ^h+sp^)'^] 




or 


(5.11) 


1 


^ i(h+spj^) -(h+sp^ ) } r^+ {r^(hj^+sp^ ^ ~'^2 , 

s= ^ \ ^ ^ -* 


^^ 1-^2 ^ 


(5.11) 
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S.2>l.Axial Force I 


The axial force is given by 

% ^2 

^3 = 2 // ® 


(5.12 } 


On sxibstitution of the value of p from equation (5.11) into 

equation (5.12 ) we obtain 

^ rg C{(h+spj^)--(h+sE^)^!r^ . 

= V / ^ *• — 

^ 1 


^^1-^2 ^ 


r^(h+sp^ -h(rjj^-r2 ) 


(ri-r2 ) 


1 r d8 dr 

m 


/ [n-<h+sp 2 _)^-(h+sp^ )^ } r 


m 


- r2^(h+sp^ ^*^”^2 }}^'^^-(r2 )h] 


de 


Now 


h dQ =/””■ (Ctr r^^ cos0)d0 = cn 
0 o 


(5* 13) 
(5.14) 


Using equation (5.14) in equation (5.13)/ we get more simplifed 


form for axial force as I 


2r 


F =- c, 

z s 


m 


71 4 4 

[ 'n;c(r-, - r, ) -/ (h+sp, ) -(h+sp^ ) } 

X z 0 

+ { r2^(htsp2 )^-t2 (h+sp^)”^ H ]<i6 ] (5.15) 


on further simplification F^, becomes 


Fz = 


2rjj^(r^ -r^ ) 


[ - Ttc t r J (htSp^)V(htSp2 )'^ 1 d0 ] 

^ * ° (5.16) 
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5*2 .1.1 Axial Force For No Flexibilitv Case J 

The axial force for the case when none of the 
faces is flexible is given by 


F (s=0) =2r^(r2-r,) / 

^ "■ 2 1 ^ + h2 > 


5 »2 ,2 Radial Force • 


The radial force is given by 
n 2 

Fp, = -2/ / COS0 r de dr 


o r. 


where 


t.. = 


- - 


r 2 8r 

From equation (5.6), the value of becomes 

\4 /, . _ \4. 


c„ = 

r 


^ (h+spj^ )^-(h+sp^ ) } 

^2 


(5.17) 


(5.1'3) 


(5.19) 


8s ^ 

Using equation (5.8) and (5.9) and equation (5.18) in 
equation (5.17) and then integrating, we get 

I, cose de + r~ /% cose de (5.20) 

>s 1 3 s • .4 


■R 6s 


where 


II = (h^-hsp^)^ 


(5.21) 


and 
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ri(h2+si^)^-r2<h^+sp^) 

Jl) l | l m i— I l■I H ^I|ll ^■ lll ■ ■ ! I — n i l III . H ■— I— 

2 2 

(hj^+sp^ ) +(h2+sp2^ 


( 5.22 ) 


where 


hj = c + r cosQ 


= c + r r 2 COS0 


5.2 .2.1 Radial Force For No-Flexibility Case ^ 


When faces are not flexible, i.e. s =0 the radial 


force is given by 


F^(s=0) = 2r f 
R m ^ 


.71 

— ihp^i 


(5.2 3 ) 


5.3 NON-DIMENSIONAL SCHEME 


We resort to the undermentioned non-dimensional 


scheme 


r r, 

- . m . — 1 . 

r^ = — - / r, = -r- , e = 

m r 2 1 r 2 


— ? ; — ; F»= 

c ' R ^2^2^ 


= 1 + er^cose ; h 2 = ~ = 1+ecose ,* s = 




Using equation (5.24) in equation (5.16) and equation 
(5*19) we get the non-dimensional forms of the axial 
force and the radial force as given below ' 

2r (1-r, ),. T 71 , - .4 

f la i_[ _„+ ■ - V tt f ^ (Its r cosG +sp^ ) 

- - L TIT 1 / 4 rn ■"! 

s ^ 

( 1 + e cose + i>^ die] 


(5.2 4 } 


(5;2 5 ), 
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and 

1 r'^ - 1 - 

Fr = "*z — J 1 . cose de+ I- cose de (5.26 

^ 6s 0 ^ 3s O 2 

where 

~ (h2 +s ) ^ (5 *2 7 

and 

- ( 5 - 28 } 

(h^-f«pj^ -hii^s 

Further, we calculate non-dimensional axial and radial 
forces abinitio to find expressions for these forces for a 
case where no elastic deformation occurs, i,e. when s = 0. 
Thus 


Pg (s=0) 


5 — = 2r (1-r. ) / t-^ — :: ) 


de 


(5.2 9) 


and 


= , rR(s=o) n 

F_ ( s =0 ) =£ - 2 r J de 

^ - o (^f^) 


(5.30) 


Now we define the following quantities I 


P -F (s=0) 

F = { ■ — 5 — J } xlOO 

^ F (s=0) 

z 


(5 .31 


and 
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F -F (s=0) 

Fp = { } X 100 (5.32 ) 

Fj^(s=0) 

From definitions in equations (S29) and (5.30) 
it is clear that F^ and 1^2 refer to the percentage difference 
in the values of the axial and the radial forces respectively^ 
in the present case from their corresponding values when 
s = 0. 


5.4 RESULT AND DISCUSSION 


The values of F^ and F 2 are obtained by integrating 
equations (5.31) and (5 .32 ). numerically by Romberg method. 

Both Fj^ and F 2 are functions of p^ = (inner to outer 

pressure ratio)/ r, = r— (radii ratio)# s (elastic parameter) 
and e (tilt parameter). The curves for F^^ and F 2 are 
plotted against the tilt parameter e (figures 5.1 to Si4:)^ 

The value of r^ is fixed at 0.9. 


Figure 5.1 represents the graph between F^ and e 
at a fixed value of p^, i.e., p^ = 0. This figure indicates 
that the percentage difference increases slowly with 
an increase in the value of tilt parameter when it is 
less than .8 (approx, ) but for e > .8# the rate at which 
F^ increases with e is faster; The three different curves in 
figure 5.1 make it clear that the value of axial force 
increases with an increase in the value of elastic parameter 
s. Since s = # an increase in s may either be due to a 
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decrease in the value of the clearance c of the centres 
of seal faces or an increase in the value of s or both. 
Initial decrease in the value of c will ccane about 
following the shifting of the flexibly mounted 
face towards the secondary seal ring whereas an increase 
in the value of s is done by a more flexible nature 
of the face material. This happens due to a decrease 
in the modulus of elasticity of face material. It can 
thus be interpreted from figure 5.1 that for a seal 
with a more flexible face material , the percentage 
difference in the value of axial force is higher. 

Figure 5.2 corresponds to the case when p^ = lOOO. 
This indicates that the pressure at the inner boundary, p^ 
is higher than the pressure at the outer boundary, P 2 * 
Whereas in figure 5.1 is less than • This 

figure 5.2 depicts a behaviour of similar to that 
of Figure 5.1, i.e., the percentage difference in the 
axial force increases with s. However, there are 
significant quantitative differences. For pj_ = O, the 
maximum percentage difference is of the order 14 o/o, 
whereas for pj^ = 1000, this difference goes as high as 
60 o/o , Also, for = 1000, the behaviotrc of F^ 
with e is quite different from that of F^ for p^^ = 0 
F^, for pj^ = 1000, remains almost constant with e till 
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e < ,6 (appicK. ), But decreases with a further 

increase in e. 

Each of the figures 5.3 and 5.4 represents curves 
of the percentage difference in the radial force, , 
versus the tilt parameter e, for various values of s. 

Figure 5.3 is a graph with p^ = 0 . This figure reveals 
that the value of F 2 is less than zero and they show a 
very little variation v/ith s. This is a very remarkable 
result in that it indicates that the flexibility of face 
material does not have a very significant impact on the 
radial forces. Moreover , a negative value of F 2 indicates 
that the radial force infact decreases when seal faces 
become flexible. Thus in view of [ 77,78 ] the flexibility 
reduces slightly the inward pumping that is generated by 
the radial force. Also, when e increases the radial 
force increases (from fig. 5. 4 ) which is in perfect 
agreement with the behaviour of radial force with 
misalignment shown by Etsion and Sharoni [ 77 j and Etsion 
[ 78 ]. 

Figure 5.4 is a graph between F 2 and e for p^^slOOO. 
The qualitative behaviour is reversed as compared to that 
in figure 5.4. This figure indicates that F^ increases 
as s increases which implies that the radial force increases 
with the flexibility of the face material. However, the 
percentage difference is not very high. Also, for every 
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value of s, the percentage difference decreases 
with an increase in the angular misalignment parameter e . 

The hydrostatic component of the radial force 
acts along the line BB which connects the highest 
and the lowest points of the primary seal ring. When 
the high pressure is on the inside periphery of the 
seal (as is the case presented in the figixre 5.4), the 
radial force F^^ is directed towards the point of maximum 
film thickness. On the other hand, when the high 
pressure is on the outer periphery, the radial force 
Fj^ acts in the direction of minimxom film thickness 
(this has been presented in the figure 5.3). Thus 
the percentage difference ^2 opposite behaviour in 

both the cases (fig. 5.3 and 5.4). 

Radial force ’ effects appear to be more 
pronounced when the high pressure is on the outside 
periphery of the seal [78 This is because such an 
arrangement is inherently unstable [ 65 ] due to 
nonrestoring hydrostatic moment which tends to increase 
any angular misalignment and thereby raise the magnitude 
of radial forces. From figure 5*4, it is clear that 
in case of higher pressure on the outer periphery radial 
forces are decreased because of flexibility of a. seal 
which implies that the damage to the stability of the 
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seal is lessened* Moreover the radial displacement 
of the primary seal ring becomes less following which 
the leakage through the secondary seal, which is generated 
due to inward pimping in such cases, is also lowered* 

It may thus be concluded that rendering the 
faces flexible is in the interest of stability of a 
radial face seal which is having high pressure on the 
outer periphery and operating under indequate separating 
force. 



S9 



Fig. 5.1 Variation of the percentoge difference in the axiai forces 
with the flit porometer for various volues of elasticity 

parameter 
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CHAPTER VI 

HYDROSTATIC ERESSURB EFFECTS IN A MISALIGNED RADIAL FACE 
GS^L WITH A HON-NBWTONIAH FLUID 

6-1 INTRODUCTION 

Akial forces in excess of those theoretically 
predicted for aligned flat seal faces are now well 
established [ 8,27,42 ]. These axial forces are 
crucial in providing the mechanism that separates 
the seal faces and controls the gap between thm. 

Several different mechanisms have been proposed to 
understand the origin of these forces. Of these 
mechanisms, hydrostatic pressure ''hich is most 
important- 

Several investigators, including Cheng and Snapp [23 ]# 
Snapp and Sasdelli [ 43 1# Metcalfe L 44 1# Mayer L 45 1# 
Metcalfel 56 ], Etsion [ 64,65 ], Lebeck [ 70 ], Etsion 
[73 ], Etsion and Burton [ 75 ] and Etsion L 76 ] have 
refined the techniques of determining the radial face 
profile and the hydrostatic pressure distribution. In 
most of the earlier investigations, the lubricant was 
treated essentially to be a Newtonian- However, the 
non-Newtonian behaviour is almost invariably observed 
in most lubricating processes. This may be n consequence 
of the fact that the fluids such as plastics, greases. 
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slurries etc. are being produced indurtrially in 
increasing quantities and these fluids are finding 
use in the liibrication systems. This non— Newtonian 
behaviour of the l-ubricant may also become prominent 
because of the contamination and the other severe 
operational requirements. In such situations, the 
Np>wtonian postulate, which assumes a linear relationship 
between the shear stress and the rate of strain, is 
violated. 

Various theories have been developed in recent 
years to describe the flow behaviour of a non-Newtonian 
fluid based upon the ftindamental constitutive relation 
that the shear stress is equal to the product of 
consistency and the rate of strain where the consistency 
itself may be a function of both shear stress and rate 
of strain. One such model is that of the 'power law' 
fluid in which shear stress varies as some power of 
strain rate. The power index, say n, is of immense 
significance in the study of the power law fluid flow 
problons. Its deviation from vinity determines the 
change in the Newtonian behaviour of the lubricant. 

For n = 1, the fluid is Newtonian ,* for n < 1 , it 

characterizes a pseudoplastic behaviour and for n > 1, 

it is a dilatant fluid. 
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The value of the flow behaviour index and consitency 
for some non -Newtonian fluids which obey the power law model 
listed below in table 6*1 Ref.[ 88,89 

Table 6.1 Typical values of the flow behaviour index and 

■ consistency variable for some well known non— Newtonian 

fluids » 


Composition of fluid Flow behaviour Consistency 

index 


23.3 per cent Illinois yellow clay 

0.22 9 

0.863 

in water 



67 per cent Carboxy-methyl -cellulose 

0.716 

0.121 

( CMC ) in water 

- 


3 . 0 per cent CMC in water 

0.566 

2 .80 

23 per cent lime in water 

0.178 

1.04 

4 per cent paper pulp in water 

0.575 

6.13 

54.3 per cent cement rock in water 

0.153 

0.331 

14.3 per cent clay in water 

0.350 

0.0344 

25.0 per cent clay in water 

0.185 

0,204 

40.4 per cent clay in water 

0.132 

2.30 

27 per cent Titanixm dioxide water 

1-41 

Data not 

suspension (0.2 -1.0 micron, particle 

size) 

available 
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In recent years^ this non-Newtonian model has 
found extensive application in various lubrication 
problems. See reference [ 88 ] for detail. However, 
it has not been applied to the case of a misaligned radial 
face seal. Tanner [ 10,11 ] did study the effects of 
non-Newtonian fluid for seals but his model dealt with 
the viscoelastic fluids without considering the misalign- 
ment of the faces . 

It is thus likely that the study of radial face 
seals when coupled with the concepts of misalignment 
and flexibility of faces may lead to significant non- 
Newtonian effects. Hence the study of these non -Newtonian 
effects through the power law model may yield additional 
information which will help in the designing of face 
seals. 

In this chapter, therefore, power law model is 
used to investigate the effects of misalignment andflow 
behaviour index on axial forces and tilting moments. 

The curvature effects that were neglected by Etsion [ 65 ] 
are also given due considerations. Seal characteristics 
such as axial force, tilting moment and side leakage 
are obtained as functions of the flow behaviour index 
and the tilt parameter.. Apaart frcm this, a comparision 
has also been performed between the Newtonian and non- 


Newtonian cases 
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6.2 MATHEMATICAL ANALYSIS 

In this section, the mathematical analysis for 
hydrostatic face seals is developed considering the 
lubricant as power law fluid. The schanatic of the 
misaligned face seal is given in figure 3-1. It is 
assumed that the power law fluid enteres between the 
faces, whose inner radii is rj_ each, at a pressure p^^. 
h is the film thickness which depends upon the angle of 
tilt y and the variables r and 0 . 

We shall first derive the generalized Reynolds 
equation applicable to the model under consideration 
and then apply this to study the seal behaviovur. 

The stress strain relation for the motion of a 
power law fluid is given by 

Cij=m ; i,3=l,2,3 (6.1) 

(The origin of coordinates system lies at c/2 ) 
where c . . is the stress tensor, e . • is the rate of 
deformation tensor, m is the consistency and n is the 
flow behaviour index. The consistency m, in general, 
can be a function of both shear stress and shear rate. 
The consistency may be interpreted as a parameter 
depending upon the viscosity of the fluid. For a 
Newtonian fluid ( n=l), the consistency is identical to 
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the fluid viscosity. The dimensions and the numerical 
values of m depends upon n. 


to 


For one dimensional case, equation 


(6.1) reduces 


, dv- n-1 

The equation of motion and the equation 
ing the flow of incompressible fluid, 
case, one given by 


( 6.2 ) 

of continuity govern - 
in one-dimensional 


ep 

3F— 


= U 


(6.3 ) 


i- = o 


( 6.4 ) 


and 


1 

r 




(6.5 ) 


respectively. 

Substituting for c from equation (6.2 ) in equation 
(6.3 ), we get 



Since the flow that occurs is pressure induced, 
with a constant pressure in the axial direction, the 
flow will be symmetrical about the film thickness. In 
view of this the boundary conditions can be stated as 


follows 
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z 


= 0, 



0 , 


z 5= h/2^ V = 0 and w == 0 


(6.7) 

( 6 . 8 ) 


Integration of equation (6.6), using the boundary 
conditions (6.7) and (6.8), yields the velocity distribution 
V, in the small gap between the mating faces of seal, as 
follows 


= ( 


n+1^^ 


1 9P )l/nr 
m 9^ ^ 


fSii) / n+l X 


(6.9 ) 


Substituting for v from equation (6.9) in equation of 

continuity (6.5) and then integrating we get 

1 ^ 

i J . i J£L) ^ (-) ^1-0 

r ar 2n-fl^^ m dr 2^ j - u 

For isoviscous fluids, equation (6.10) reduces to 

1 


/2 n +1 \ - 

-^L r h ^ dr ^ -1 = ° 


Integrating equation (6.11) we get 


( 6 . 10 ) 


( 6 . 11 ) 


p = B(e ) / 


dr 


r, n , (2n+l) 
1 r n 


•f- C (6 ) 


( 6.12 ) 


where B and C are the constant of integration and are 
probably the functions of 0 only. 

The film thickness h for a misaligned radial face 
seal is 


h = c t 7 r cos 9 


(6.13 ) 
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The pressxare boundary conditions as per the situation 
quoted earlier are 


p = p^ at r = ( 6 . 14 ) 

p = O at r = r 2 (6.15 ) 

Substituting for h from equation (6.13) ib equation 

(6.12) and integrating, using the boundary conditions (6.14) 

and (6.15), we get 


= B(e) / 
^1 


dr 


r^(c + y r cos6 ) 


■('fn+I'J (6.16) 


where 


C(0 ) = p^ 


(6.17 ) 


and 


B(e) = 


~ Pi 


dr 


r2 _ 

^1 r^(c+r r cos 6) ^n+l) 


For n = 1 


/ : 

r(c+')' r cose)' 


p = -PjL 




dr 


r(c+‘)' r cos8 )' 
6.2.1. Axial Force * 


(6.18) 


(6.19) 


The axial force, v/hich provides separation 

mechanism, is given by 
2% r 

F^ = / / p r de dr 

^ o ^1 


( 6.2 0 ) 
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Integrating this equation by parts, we obtain simpler 


formula for F_ as follows 
z 

- '” 4 % - r / / ^ ^ '»® 

w X. ^ 

Supstituting for from equation (6.12 ), we get 


( 6 . 21 ) 


X 


■2 dr 


- 7t 


2 1 r -"1 


2% •b, „(n-2 ) 


(c+_2lrcqs0) 


(2 n+1 ) 


dr 


/ 

r^(c+r r cos6 


( 6.22 } 


For Newtonian case (n=l) 


r ^ r dr 

2tc -/ r 2^ Tct7~r~cos9~7^ 

F^(n=l) = - n <— J ) (6-23) 


/ 


dr 


o 

(c+'y r COS0) 


6 *2 .2 Tilting Moment I 

Tilting moment, the moment about x axis, is given 
by (also called the restoring moment) 


2n Tg 2 

M p r cose de dr 

^ o 


(6.2 4 ) 


where (-) sign is used in order to obtain the positive 
value of restoring moment. 

Integrating equation (6.24) by parts we get 


1 .2-n: ap 3 

M.. = y r cose dr d6 

o 


X 


(6.25 ) 



102 


dp 

Substituting for from equation (6.12 ) in equation 


C 6 .2 5 ) we get 


dr 


^ ( n-3 ) ^ c+r r cose ) ^ 

/ (“T ) cose de 


dr 


^ o" j'^2 

^1 r^(c + y r cose 


(6.26) 


e *2 «3 * 


Q 


The leakage across the seal is obtained from 

^ 2 n+1 j 

- 2 "" (i^)Vn^^e (6.27) 

J 2n-H 9. m fir 


Substituting for — from equation (6.12 ) in equation 
(6«27 ), we get 


f n ^ , ^1 ^ 1 /n 
^ ^ (2n+l]^ 2m ^ q*^ 


de 


(/ 


dr 


a/n 


^1 r^Cc+y r cose ) ^ 


(6.28) 


6.3 NON-DIMENSIONAL SCHEME 


For nondimensionalizing the various quantities 
involved, the undermentioned scheme is adhered to 


y r 




1 • !* 


"2 . - 


^ r2 = — . r 




^1 ' ^ ^2 =2 


F 3M 

— Z . P • _ k 

F = — * — , p = — ^ , M_. = 


X „ 

Pl^l 


(6.2 9) 


Thus the dimensionless pressure. 


force and tilting 
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moment assume the following forms 


/ 


dr 


P 


^ r^(l+ e r cos 6) 


(2 ii+l ) 


( 6 . 30 ) 




dr 


“U / -», ** ^ nH"! } 

r (1+e r cos B; 


r^ dr 


f = --Ft + p / C 


r a t 

„ In -n / , ^ ^£,>(2 n+1 J 

2it 1 r (It e r cosG; 


“) d.9 { 6 • 3 1 J 


dr 


/ 

^ ?^(l+e r cose)^^"^*^^ 


and 


M 


X 


r ^ 

a-rt 1 “(n-3}Q^g. - ^ 

J* ( ... . — — )cos0 d6 (6.32 ) 


dr 


/ 


(2 n+1) 


^ r^(l+e r cos6 ) 

where e = 0, i-e. for parallel seals 


{ e =0) = -Tt + ^ ^ 

M (e = O) = 0 
X 

Now when e = O and n = 1/ i-e. for parallel seals 
Newtonian lubricant 

n -3 ) 

F(n=l,e=0) = 

2 log r 2 


(6.33 ) 

( 6 . 34 ) 

with 


(6.35 ) 
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Nondimensional scheme for the leakage is complicated 
for power law model because the dimension of the consistency 
m depends on the flow behaviour index n. Therefore# 
we consider the ratio of tilted faces leakage to 
corresponding parallel faces leakage to study the 
behaviour of leakage. Thus we denote the ratio by 


Q = 


Q 

Q(e=0) 


% 


f 

o 


de 


( / 


dr 


r’^Clta r cose)^^^"^^^ 


)Vn 


. . r 

o ( 

^ 1 


de 


( 6.3 6 ) 


dr - 




n 


6.4 RESULTS AND DISCUSSION 


Equations (6 .3 1 )and (6.32) give expressions for 
the dimensionless axial separating force F and the 
dimensionless tilting moment M respectively. Each 
of them is a function of the flow behaviour index n 
and tilt parameter e. Equation (6.33 ) represents the 
expression for separating axial force when e = 0 and 
is positive. Equation (6.34) depicts, as expected, that the 
tilting moment vanishes for e =0. 

These equations have been solved ntmerically using 
Romberg method, for inner-to-outer radii ratio equal 
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to 0.9. This is in accordance with Etsion [ 65 ] in that 
the author mentions that for all practical purposes the 
radii ratio is greater than 0.8 . The other two important 
parameters are the flov/ behaviour index n which is made 
to vary from 0.2 to 2 .0 [ 89 ] and the tilt parameter 
which is made to vary from O to 1. Each of the quantities 
has been plotted against both the parameters n and e « 

Figure 6-1 is a graph between the non-dimensional 
axial force F^ and the tilt parameter e, for various 
values of n. The value e = O corresponds to the case 
when the faces are parallel. The figure indicates that 
as the value of e increases the value of axial force F^ 
for pseudoplastics increases negligibly till e =0.6 (approx.). 
When £ 0.6 F_ increases more rapidly. For dilatant fluids^ 

F^ decreases first reaches a minimum aroxond e =0,7 and 
then increases . 

Figure 6.2 shows the variation in the axial force 

F with n. This figure shows that for all values of e 

the value of the axial force F decreases as the flow 

z 

behaviour index n# increases. 

F 

•^z 

We define a ratio » i.e. the ratio of the 

F^(£=0) 

axial force for misaligned radial face seal to the axial 
force for parallel face seal. This ratio has been plotted 
in figure 6.3 and 6.4 against e and n respectively. It is 
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seen that for dilatant fluids (n>1.6 approx.) this ratio decr- 
eases with e first reaches a minimum around e = 0,7 and 
then rises sharply. This indicates that a small angle of tilt 
has a determimental effect on the axial force. For fluids 
for which n < 1.6, the ratio is always greater than unity 
which signifies that the separating force is always larger for 
a misaligned face seal than for a parallel face seal. 

Similar interpretations can be made from figure 6.4. 

Prom these figures it can be inferred that for a face 
seal that \andergoes massive wear because of insufficient 
generation of separating force the use of pseudoplastics 
becomes favourable in that it increases these forces and helps 
in maintaining the faces apart. But for a face seal that 
starver for lubrication because of frequent huge leakage, 
the use of dilatant fluid will be beneficial. 

Figures 6.5 and 6.6 represent the graphs of the tilting 
moment M against e and n respectively. These figures 
together reveal that the tilting moment increases as 

the value of e and the flow behavio\ir index n increase . 

For a pseudoplastic fluid, is much less than that for 
a Newtonian and a dilatant fluid. Thus the dilatant 
fluids are more useful where face seal stability is more 
frequently hampered. Because these fluids produce more 
restoring effects on the primary seal ring. 
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Figure 6.7 is a graph of the leakage ratio defined 
in equation (6.36) versus the tilt parameter e. This 
ratio increases as n decreases. The rate of decrease 
is much more rapid for small values of n and large values 
of e# i.e. for pseudoplastic fluids with large angle of 
tilt. This suggests that the leakage is larger for pseudo- 
plastic fluids than for Newtonian and dilatant fluids. 

This may be the consequence of increase in separating 
force and decrease in restoring moment for pseudoplastics 
as indicated in figures 6.1 and 6.2 . 

On the basis of above findings it may be concluded 
that for a radial face seal, the dilatant fluids as the 
lubricant are more useful than the Newtonian and pseudoplas- 


tic fluids 
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R.ADIA.L AND FRICTIONAL FORCES IN A MIS7vLIGNSD RADIAL FACE SEAL 

WITH A MON-NEWTONL'II FLUID 

7 . 1 INTRODUCTION 

Ever since Denny [ 8 ] demonstrated the existence of a 
lubricating fluid film between the mating faces of a mechanical 
face seal, a considerable amount of fundamental research has been 
clone in an attempt to understand the mechanism that enables face 
separation. A recent review is by Ludwig L 54 ] in which various 
models and theories are listed. Much attention has been paid in 
the last two decades to the radial face seals. In particular 
axial forces and tilting moments were studied. The radial forces 
\ 7 hich may cause the radial displacement of the primary seal ring 
has not been analysed throughly. 

These forces may be of Important consequence in the onset 
of radial misalignment. Radial forces may well be the origin of 
radial eccentricity and hence directly influence the primary seal 
leakage. In addition to their influence on leakage^ the radial 
forces may also affect the dynamic stability of the seal. 

Recently Etsion L 78 ] has analysed the radial forces in a 
misaligned radial face seal and shown that radial forces 
are most pronounced when the high pressure is on the 

This chapter was presented in paper form at the Winter Annual 
Meeting of the -^^mencan Socrety of Ilechamcal Engxneers, Phoemx^p 
Arizona, November , 1982 . 

(To be published in Jour. Fluid. Engng», Trans . ASME ). 
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outside periphery of the seal. In such cases the inward 
pumping generated by the radial forces, adds to the 
hydrostatic leakage and this may provide a possible 
seal failure reason. 

Shapiro and Colsher L 47 ] and Ludwig and Allen [ 55 ] 
considered the dynamic behaviour of a radial face seal and 
studied the frictional force on the secondary seal ring. 

This friction provides the damping of the axial and 
angular vibrations of the primary seal ring and is dependent 
on the radial pressure exerted by the primary seal ring 
and seal housing on the flexible element of the secondary 
seal. 

Most of the work devoted to the understanding of the 
seal mechanism, which has so far appeared in the literature 
assumes that the Ixibricant behaves essentially as a 
Newtonian viscous fluid. Non-Newtonian behaviour is, 
however, almost invariably observed in various Ixxbrication 
processes. The reason for this has been amply discussed 
in the previous chapter. 

It is thus felt that the application of the power 
law model to the analyses of a misaligned radial face 
seal may well lead to a better visualization of the 
origin of the radial and the frictional forces which 
affect the seal performance considerably. The results 
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of the analysis might also be useful in an instance where a 
non-Newtonian situation is encountered. Curvature effects 
neglected by Etsion [ 78 ], might be of important consequence 
in this context and hence are considered is this chapter, 
further# in this analysis no appro jdLmat ion of the type made 
by Etsion [ 78 ] has been resorted to, thus leading to 
more accurate results. 


7 .2 MP^TNEMTiTICiL ANALYSIS 

The schematic of misaligned radial face seal is given 
in figure 3.1. 


Mathematically# for a power law fluid# 


■ 3v , n-1 6v 

The equation of motion is 

dp _ d c 
dr dz 

and the equation of continuity is 


(7.1) 


(7.2 ) 


1 d / X ^ 

r ar az 


= 0 


(7.3 ) 


Substituting the value of v from equation (7.1) into 
equation (7 ,2 )# we get 


dr 





Ne confine our attention to the region where 
Thus equation (7.4) reduces to 


(7.4) 



< 0 . 


Integrating 


equation (7.5) with the boundary conditions 
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dv , 

^ — = 0 at z 


Yi/l 


(7.6 ) 


V =r O at z = h and 2 = 0 


VJe obtain the velocity profile as I 

n+1 n+1 


,h, n 

V = ( [ -2^ 


•(z-h/2 ) 


n 


m 


ntl 

n 


1, 2 <_ h 
(7.8) 


and similarly 

n-fl Itn 

1 do ' r ^ "" -(h,>2-z) n 

= ^ i df”^ 3' Oi 2 i 


ntl 

n 


(7.9) 


where h is the film thickness. 

Substituting for v from equation (7.9) in the 
equation of continuity (7,3) and integrating with respect 
to z with the conditions 

w = 0 at 2 = h and z =s 0 (7.10) 


we get the final equation controlling the motion of a 
power law lubricant in the gap between the faces of a 
radial face seal as below I 


2 n+1 

i 2n /hx n . 1 

r dr ^ 2 n+1 ^2 ^ m dr 


1 

n 


= o 


(7.11) 


The boundary conditions for the pressure are I 
p = p^ at r = rj_ > P - P 2 ^ = ^2 


(7.12 ) 
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where is the inner radius and r 2 is the outer radius. 

Integrating equation (4.11) twice and then using 
the boundary conditions in the equation (4.12 } and 
assiaming the consistency to be constant, we get I 


ar 


r h 


Pi 


where 


0 ,( 0 ) = — - 
1 r. 




The film thickness h is (see figure 4.1) 

h = c + r r cos© 


(7.13 ) 


(7.14) 


(7.15) 


where c is the axial clearance and r is the angle of tilt. 
Thus p has the following form \ 


[- 


<E^-Pl) 


-] X [/ 


dr 


dr 


- 

'^1 r^(c+ 7 r cos 


r^(c+r r cos©)^^'^^ 


(7.16) 


3+Pi 


(i) Radial Force '* 


The radial force in a seal is defined as 
% 


F = - 2 / / 

o 


t: cos© r d© dr 


(7.17 } 


Thus the radial force F^ on the surface z = h is obtained 
by substituting the value of c at z = h, in equation (7.17 ). 



220 


Thus 


% 




dr 


/ — 

^ r^Cc+r r cos 9) 


'? n+1 


{ / 


dr 


^1 r^'’^(c+‘y r cosG)^^ 


} ] co^ d0 (7.18) 


(ii) Frictional Force I 


The radial frictional force is defined by 
TX ^2 

P = - 2 . f f crd0dr 
R o, Ti 


(7.19 ) 


and at 2 = h the radial frictional force assumes the form 

% 

I 

o 


r 

Fr = (E^-p^)cr2 / LC- 




/ 


dr 


r^(c+y r cose)^^'*’^ 


X { / 


dr 


1 r^~^(c+Y r cos0)^^ 


} ] d9 (7.70) 


7.3 NON-DIMENSIONAL SCHEME 

We resort to the following dimensionless scheme .1 

y r„ r, 

- .- r*^ 2 * - .1 

n = 1 + e r cos© / r = -p- / £ = ' ^1 ~ r 

2 2 


F 


^ = 




(7.21) 


"^1 ^^^^2 -Pi ^*^^2 

Thus the non-dimensional radial force becomes 
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= i I cose de 

o 


(7.22 } 


wnere 


I = L 


/ 


dr 


3x[/ 


dr 


r^d-i-e r cosO)^^"^^ 


r^^^Cl+e r cose)^^ 
(7.23 ) 


And the non-dimensional frictional forces is 


Tl 


= / L X 
° / 


dr 




dr 


r^(l+e r cosB)^’^'^^ 


— -n-1 , T „ \2 n 

r (1+e r cose; 

(7.2 4) 


]d6 


For Newtonian, fluids, i.e. when n = 1 


F^(n=l) = 


1 TC 'll ( 9 ) 

-y ^ ae , . 0 

O 2 


(7.25) 


where 


( 0 ) = [ 


1+e cos e 


1+e cose 


(7.26) 


and 

q2 ^ 


1+e r-, COS0 

[log(^;; ^ ) + ( 


r, +c £-tCooO 1+e cos 6 1+ e r, cos0 

11 


. If 1 


(7.27 ) 


(1+e cose) (1+e r^ cose) 


*^nd 


n q^(0 ) 


FR(n=l) (0 ) 


(7.2 0) 
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For n = 1, e = 0 , 

F^(n = 1, e = 0) = 0 (7.2 9) 

il-r, J 

F (n = 1, e = O) = -Ttx ^ (7.30) 

log(rj|^ ) 

7.4 RESULTS AND DISCUSSION 

The most important parameter for discussing the 
behaviour of various seal characteristics, in the 
present context, is the flow behaviour index n. The 
fact that the power law model can actually be used for 
realistic fluids is obvious from table 6.1. It is seen 
that the value of n, generally lies between 0.5 to 2.5. 

In the discussion that follows each of the seal character- 
istics is analysed from this view point. 

Fig. 7.1 is the graph for non-dimensional radial 
force versus the tilt parameter e . It is seen that 
an increasing value of the tilt parameter increases 
this force. This force is almost a linear function of 
e for all values of radii ratio ^ 2.^2 flow behaviour 

index n. The force F^ acts along the line connecting 
the maxim\jm and the minimxm film points of the primary 
seal ring. The direction of this force is determined 
by the sign of high pressure is inside 

the periphery of the seal, directed towards the 

point of maximxm film thickness, and vice versa. It is 
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also seen that lower the value of n, higher the radial 
force . 

Fig, 7.2 is a plot of friction force Fj^ vs, tilt 
parameter. It is seen that higher the value of tilt, 
lov7er the value of frictional force on the primary seal 
ring. Thus tilting may be utilized in the reduction 
of the frictional forces. However, the tilting does not 
produce a significant difference in the frictional force. 
The behaviour of the frictional force, with respect to n, 
is similar to that of the radial force. 

The qualitative behaviour of Fig. 7.3 is similar 
to that of Fig, 7.2 with difference that Fig. 7.3 gives 
the value of F^ as a function of n, thus includes a wide 
range of fluids for different tilt parameters and radii 
ratio rj^ / ^2 ‘ curves for n < 1, depict an increase 

in the radial force (Fig. 7.1) , frictional force (Fig. 7.3 ) 
It can thus be concluded that for pseudoplastic fluids 
the radial force as well as the frictional force increases 
when compared with the corresponding Newtonian results 
(n = 1) . 

Table 7.1, given a comparision between the 

radial forces M (obtained by Etsion L 78 ], neglecting 
o 

curvature effects) and Mj^(present analysis, curvature 
effects included) for various values of eccentricity ratio 
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and inner to outer radii ratio, r^, for Newtonian fluids. 
It is seen that the curvature effects are not significant 
for values of r^^ exceeding 0.9. However, for lower values, 
the inclusion of curvature may decrease the radial force 
by over 2 percent in the Newtonian case. A similar result 
may be expected for non-Newtonian fluid. 
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CHA.PTER VIII* 

EIASTOHYDROSTATIC LUBRICATION OF CIRCULAR PIATE THRUST 
BEARING WITH POWER LAW LUBRICANTS 


0.1 INTRODUCTION 

The increasing niimber of practical bearing employing 
layered elastic solids as mating surfaces has led to a 
substantial interest in the elastohydrodynamic and elasto- 
hydrostatic lubrication of such bearings. To avoid wear and 
seizure and to make the supply of the liibricant easier, 
thin layers of softer materials, such as the stiffer plastics, 
are sometimes coated on the bearing surfaces. In such cases 
the bearing characteristics and film geometry are modified 
substantially in comparision to the corresponding rigid 
bearing. 

Dowson and Taylor L 90 ] analysed the ' behaviour of 
elastohydrostatic circular thrust plates. An increase in 
load carrying capacity and change in film geometry which 
was predicted theoretically was verified ej^erimentally in 
the same paper. Although the theoretical analysis was based 


^This chapter has been published in the form of following 
research paper ''Elastohydrostatic Lubrication of Circular 
Plate Thrust Bearing With Power Law Lubricants", Jour.L\ib. 
Tech., Trans, ASliE, Vol.104,19^ ,p.2 43 « 
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>» □ 
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schematic diagram of elastohydrostatic thrust bearing 
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upon some approximations, some useful deductions were made 
for such bearings. 

The results of [ 90 ] are based on the fact that the 
lubricant used was Newtonian. This chapter extends this 
work to lubricants which may display non-Nev/tonian behaviour 
due to the presence of certain additives. The non-Newtonian 
model chosen here is that of a pov/er law which has been 
considerably discussed in Chapter 6. 

The circular plate thrust bearings considered in 
this chapter are of the type used in the fluid sealing 
end are known as radial face seals. The analysis and 
the results discussed in this chapter may therefore be 
extended to predict the behaviour of radial face seals . 

G,2 MA.THEMATICA.L ANALYSIS 

The elastohydrostatic lubrication of circular thrust 
IDlates is analysed using a non-Newtonian power law lubricant. 

Consider the flow of a power law fluid between two 
externally pressurized circular plates separated by a 
constant distance h apart. The liibricant is supplied 
through a central hole of radius r-j^ with constant pressure 

The upper and lower plates are lined with material 
much softer than the rigid backing ond the thickness of 
each layer is f/ 2 where f « r 2 / vrhere r 2 is the outer 
radius of the plate. Figure 8.1 shows a schematic diagram 
of such a bearing. 
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As has been derived in Chapter 6, the final equation 
governing the flow of power law lubricant in the gap between 
the plates is as following 

/2 n+1 \ 1 

| j[rH " (£)"]=0 ( 8 . 1 ) 

Let h^ be the film thickness which is constant in the 
unstressed case. For stressed case, when the pressure 
modifies the film shape, the film thickness h can be written 
as 

h' = h + sp (8.2 ) 

where sp is the local compression of the elastic layer, 
given by [ 90 ] 

g ^ _f (0.3 ) 

E' 


^vhere E is the modulus of elasticity of the soft material 
and E' is the elastic onstant and is given by 


1_ 

E" 



2v 


1-v 


( V is Poisson's ratio) 


(5.4) 


If the variation of f with radius is neglected, equation 
(0.1) can be integrated for pressure distribution and other 
bearing characteristics. 
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The pressure boundary conditions are 

p = at r » r^ (8.5 ) 

p=0atr = r2 (8.6) 

bolving equation (8.1) for pressure distribution with 
the pressure boundary conditions (G.S) and (8,6)/ we get 
for n d 1/ 


p = - [ ( 


( ! c- j_v, i(2n+2 ) , (2n42 )^ 

f(sp^+h^) ~h^ ^ ,(2n42) 

‘ TTiin + ^o 

11 ^ ""2 


(■ 


(n-l) (n-1) 

^1 ^2 


1 

(2n42 ) 


-J:.-- 1 ) -h ] (8.7) 

,(n-l) ^Cn-1) ° 


;or n = 1, 


log(-^)+(l+P 3 _-|“)^ log( -~) 


h 


P = 


[£ 


o 


} ' - 1 ] ( 8 , 8 ) 


log(-- — ) 
^1 


•2 . 1 Load Capacity 


The load carrying capacity of the bearing is 


2 2 

W = Tc r^^p^ + 2n S pr dr 


(8.9) 


for n 1»3, 
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= "4pi+ -h t { ' T-^ ^ 


, 1 1 

^ (n-l ) ^(n-1 ) 


^i^j:Tr-^: 


•) 


,/(2n+2 ) \ (2rn2 j -i 

+ h ; r dr J ~ 


% h 


,2 2 ^ 


( 8 . 10 ) 


for n = 1, 


TC h 


= 71 r^Pj^ - ~ 


,2 2 . 


logC^)+(l+i— -^)^log(— 


2 Tih _ “2 — =■ 'r 


-[/ C 


r ' 1/4 

i- r dr ] 


X *^9 

1 log(^) 

^1 


(8.11) 

Taking the limit s •* 0, the equations (8.7)- (8. 8) and (8.10)- 
(C.ll) yield the pressvire distribution and the load carrying 
capacity for rigid bearings as follows 
for n / 1, 


P = 


( 


.(n-1) _(n-l) 


-i ) ^ ^2 


■) ( 8 . 12 ) 


:or n = 1/ 


P = P, 


log(p-) 


log(rr-) 

^1 


(8.13 ) 


for n / 1#3 



135 


= 


2 2w Pi 


(■ 




1 ^ 


f 1,1 

TalHT^rni:^ " 




If 1 

2^:Tr-3T 


(n-l ) 


■) ] 


(8.14) 


for n = 1, 


n Pi 


w = 


- 2 


Cr^-4) 


log( 

^1 


for n = 3 , 


W 


- 1 ) 


»3 NON-DIMENSIONAL SCHEME 


We adopt the following non-dimensional scheme 


{8.15 ) 


o Pi Pj r ^1 ^1 _ 

= 7 tr^p^+ ^ L ( -y -l)-21og( ] (8.16) 


h = 


jl 

h 


r, = 


s = 


®Pl . 

h ' 


W 


w 


^^ 2^1 


(8.17 ) 


The dimensionless for various bearing characteristics for 
elastic layered bearings are given by 


h' = 1 + s p 


(8.18) 
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for n d 


P 


1 r >^2n42 ) , 

. i[ f zL 




■)( 


(n-l) (4-3-) 

-i ?, )+l) -1] 

-(n-l) 1 

(8.19) 


for n = 1, 


(1+s )‘^log(r ) - log(~) 

r. 


p-= Mt 

s 


logCr^ ) 


-1] 


(8.20) 


for n 1, 


-2 


V'. = r3_ - 


2 1 (s+1)^^"'^ ^-1 ) 

3/ [ (l-r^^’‘^^)+l] 


1 

(2 n+2 ) 


r dr ( 8 .2 1 ) 


for n = 1# 

.-. -2 

V; = r-, ;;;; 

S ”” 

(s-rl )^log(rj^)~log(~--) 

[ li- r (8.22) 

® log(r^ ) 

For a rigid bearing# equation (3.18) to (8,22) take the 
forms 

h'^= 1 . (8.23 ) 


for n ^ 1# 
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“(n-1 ) 

r. (n-1) 

f(n-l) ^ 

P = ( ; ; ) (8.24) 

l_j(n-l) 


for n = 


P = 


loq(r } 
log(rj^ ) 


(3.25 ) 


for n / 1,3 ^ 
0-2 

-2 


" = ^1+ 


( 1 -: 


- (n-1) 




) 


3 -n 1 


-1) - )] 

(8.26) 


for n = 1, 


W = 


(1-?^) 

-2 log(r-|^ ) 


(8.27 ) 


for n = 3, 


W 


-2 

r 


= — L (r^-1 )-2 logCr^^ )] 


(8.28) 


( 1 -?^) 


To Study the variations in load, a load coefficient 
is defined as follows 


W ( elas tohydrostatic }-77 ( rigid ) ^ j^qq 
W( rigid) 


(8.29) 


It is seen that is a function of the elastic 
parameter s and the flow behaviour index n for a given value 
of 



8.4 RESULTS AND DISCUSSION 
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The integrations involved in equations (8.2 1) and 
(0.22 ) have been evaluated numerically using the Romberg 
Method. The value of parameter r^^ has been taken as 0.05. 

-For most of the fluids the value of n lies between 0.5 
to 2.5. 

The most important parameter for discussing the 
effects of power law fluid as a lubricant on the performance 
of elastohydrostatic lubrication of circular thrust bearings 
is the flow behaviour index n. 

Figure 8.2 shows the variation of non-dimensional 
load capacity W as a function of elastic parameter s for 
various values of n. The rigid bearing case is depicted 
by the value at s = 0. It is seen that the load capacity 
increases, for all n, for a bearing with elastic layer 
as compared to rigid bearing. It is also seen that lower 
the value of n, higher is the load carrying capacity of 
bearing. 

Figure 8.3 is the graph between nondimens ional load 

capacity W(W = and the flow behaviour index n for 

various values of s . Again it is seen that the elastic 
property of the layer improves the load capacity of the 
bearing, for all values of n, whereas the effect of increasing 



139 


n reduces the rate of increase of the load capacity. 

The load coefficient which determines the 
percentage rate of increase of load, is a function of 
elastic parameter s and flow behaviour index n. Figure 
8,4 is the graph between and s for various values of 
n and figure 8.5 shows the variation of with respect 
to n for different values of s. From this graph it is 
seen that is higher for higher values of n and the 
coefficient increases with increase in s. The result 
reveals that although the load capacity is lower for 
higher values of n, when elastic layer is attached to the 
surfaces, the load capacity is increased and this increase 
is enhanced by higher values of n. 

Figure 8.6 is the graph for the nondimensional film 
profile h as a function of r. It is seen that the deviation 
of the film profile from the rigid bearing case Cs=0} is 
higher for lower values of n. This is because of the 
strong pressure generation for these values of n. The 
rigid bearing film profile is constant and is given by 
li = 1 (shown by dotted line). 

The results of L 90 ] and [ 92 ] are particular cases of 
present study and can be obtained by taking n = 1 and 
s = 0 respectively, the case of L 93 ] is obtained 
for incompressible fluids. 
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Again, it is to be noted that the increase in 
load carrying capacity obtained by [ 90 ] is more 
than that obtained in present case (n =1). This 
difference seens to be due to the approximation adopted 
by those authors. They have expanded the expression 
for the pressure gradient in the pov/ers of s(s=pf/E'h )- 
They further assume that s <<1, this enables them to 
resort to the first order approximation. Siibsequently, 
approximate analytical solution for the pressure distribution 
and the load carrying capacity are obtained. They further 
mention that the difference between the first order 
and the third order approximations in the dimensionless 
load coefficient may be 13 percent in the most 

severe conditions. However, in the present analysis, 
no such restriction , i.e. s<< 1 , is imposed. Hence 
the results of the present analysis are likely to be 
more accurate. 

The analysis of present chapter reveals that the 
load capacity increases with the elastic parameter s 
and it decreases with an increase in the value of flow 
bchvaviour index (fig- 8.2 ). When extended to parallel 
face seals, this result enables us to say that the face 
separating axial force is increased by the flexibility 
of the faces. Thus it may be inferred that coating soft 



material on seal-face improves the performance of 
those face seals which witness early failure because 
of inadequate axial force that provides separation of 
the faces , 



astatic) - 






s 







W' 




Fio S '/ 

wKd ^ » ftsm 



t j Q i Or^ ^, * 

:hc oiast;c7*v‘''nn''^®'°"‘®^' 

flow bennJini 

‘-’enaviour index. 


wiih 

Qf^Ous values of 


I 


I 





S » 0 2 
S a 0 . 

§ a 0.05 








Fig 8 4 Vor ia^ioh bv bimcrrsipAiess lood cop ffit ie 

porarncicr for VG'ious values of flow behaviouf in 




Fig 0 6 Variation of dim elision less load coefficient , 
with flow behdviq»jf inde)t for voOdus^^ 
of elastic pararoc te t / . 



soo 



147 


REFERENCES 


1. Pogg/A, (1946) 

Fluid Film Lubrication of Eferallel Thrust 
Surfaces , 

Instn. of Mech.Engrs., Vol. 155,p.l36. 

2. Brkich, A. (1950) 

Mechanical Seals Theory and Criteria for 
Their Design. 

Rrod. Engng. , Vol.2 1. 

3. Salama, M.E. (1950) 

The Effect of Macro-Roughness on the Performance 
of Parallel Thrust Bearing- 
Proc- Instn. Mech. Engrs Vo 1.163 . 

4. Wood, T.H.(1954) 

Mechanical Shaft Seals in the” Chemical Industry. 
Trans . Instn.Chsm-Engng.,Vol.32 , 

5. Zienkiewics, O.C. (1957). 

Temperature Distribution Within Lubricating 
Film Between Parallel Bearing Surfaces and 
its Effect on the Pressure Developed. 

Proceedings of the Conference on Lubrication 
and Wear. Instn. Mech. Engng. 

6. dagger, E.T. (1957). 

Rotary Shaft Seals I The Sealing Mechanism 

of Synthetic Rubber Seals Rxinning at Atmospheric 

Pressure. 

Proc. Instn. Mech. Engrs., Vol.171, 

7. Whitley, S. and William, L.G. (1959). 

The Gas Lubricated Spiral Groove Thrust 
Bearing. 

I. G. Report 28 (RD/CA ), (H.M. Stationary Office). 

0. Denny, D.F. (i960) 

Some Measurements of Fluid Pressures Between 
Plane Parallel Thrust Surfaces With special 
Reference to Radial Face Seals. 

Wear, Vol.4, p.64. 

9. Nahavandi, A. and Osterle, F.(1960). 

The Effect of Vibration on the Load Carrying 
Capacity of Parallel Surface Thrust Bearing. 

ASME Paper No.60-Lub-3. 



148 


10. Tanner, R. I. (1960) 

Non-Newtonian Flow and the Oil Seal 
Rroblem. 

J.Mech.Engg.Sc.,Vol.2, p.25. 

11. Tanner, R. I. (1961 ) 

The Riener Centripetal Effect in Face deals. 
Rroc.lntl. Conf. on Fluid Sealing, BHRA., 
cranfield, U.K.,MK43, OAJ, Paper El. 

12. Davies, M.G. (1961) 

The Generation of Lift by Surface Roughness 
in a Radial Pace Seal. 

Inti. Conf. on Fluid Sealing, BHRA, Paper C4. 

13. Ishiwata, H. and Hirabayashi, K. (1961) 

Friction and Sealing Characteristics of 
Mechanical Paco Seals. 

Inti. Conf. on Fluid Sealing, BHRA, Paper D5. 

14. Huhn, D. (1961) 

Theory of Fluid Sealing. 

Inti. Conf. on Fluid Sealing BHRA, Paper D2 . 

15. Mayer, E. (1961) 

Unbalanced Mechanical Seals for Liquids. 

Inti- Conf. on Fluid Sealing, BHRA, Paper E2 . 

16. Mayer, E. (1961) 

Leakage and Friction of Mechanical Seals 
With Special Consideration of Hydrodynamic 
Mechanical Seals. 

Inti. Conf. on Fluid Sealing, BHRA, Paper E3 - 

17. Summers-Smith, D.(196l) 

Laboratory Investigation of The Performance 
of a Radial Face Seal. 

Inti. Conf. on Fluid Sealing, BHRA, Paper Dl. 

18. Whiteley, S, and William, L.G. (1962 ) 

Principles of Gas Lxibricated Shaft Seals. 

Jour, of Mech-Engng.Sc. Vol.4, No. 2, 

19. Davies, M.G. (1963 ) 

The Generation of Pressure Between Rough, 

Fluid Lubricated, Moving, Deformable Surfaces. 
Lub.Engng. , Vol.2 4,No*ll . 



149 


2 0. Lyman, P.A. and Saibel, E. (1963 ) 

Leakage Through Rotary Shaft Seals. 

ProG. Fourth Applied Mechanics Congress, 

Paper Cl, 

21. Nau, B.S.(1967) 

Hydrodynamic lubrication in Face Seals. 

Third Inti. Conf. on Fluid Sealing, BHiE^A, 

Paper E5 . 

22. Kojabshian, C. and Richardson, II. FI. (1967 ) 

A Micropad Model for the Performance of 
Carbon Seals. 

Third Inti. Conf. on Fluid Sealing, BHRA, 

Paper E4 . 

23. Cheng, H.S. and Snapp, R.B. (1967 } 

A study of Radial Film and Pressure Distribution 
of High Pressure Face Seals. 

Presented at the Third Inti. Conf- on Fluid 
Sealing, Cambridge, BKRA, Paper No.E3. 

24. Findlay, J.A, (1968) 

Cavitation in Mechanical Face Seals. 

J -Lub. Tech- , Trans . ASME, p.356. 

25. Anno, J.N. , Walowit, O'.A. and Allen, C.M,(1968) 

Microasperity Lubrication. 

J. Lub. Tech. , Trans . ASME, Vol . 90, p.35l. 

2 6. Cheng, H.S ., Chow, C.Y, and Wilcock,D.F. (1968 ) 

Behaviour of Hydrostatic and Hydrodynamic 
Non Contacting Face Seals. 

J . Lub . Tech . , Trans .ASHE , Vol . 90, No .2 , p . 5 10 . 

27. Pape, J.G.(1968) 

Fundamental Research on a Radial Face 
seal. 

Trans. T-iSLE, Vol. 11, Mo. 3 , p.302 . 

28. Sneck, H.J.(1969) 

The Misaligned, Eccentric Face Seal. 

J-Lub.Tech, Trans.ASMS, Vol. 91, No. 4, p. 695 . 



150 


29. Sneck,H.J. (1969) 

The Eccentric Face Seal With Tangentially 
Varying Film Thickness. 

Presented at the Fourth International 
Conference on Fluid Sealing held in 
Conjunction with the 2 4th ASLE iAnnual Meeting 
in Philadelphia, FICPS Preprint No-lSA. 

30. Anno, J.N. , Walowit, J-A. and Allen, C.M.(1969} 

Load Support and Leakage and Microasperity- 
Lubricated Face Seals. 

Presented at the Fourth International 
Conference on Fluid Sealing held in 
Conjunction With the 24th ASLE Annual Meeting 
in Philadelphia, FICFS, Preprint No. 21. 

3 1 . Findlay, J .A , ( 1969 ) 

Measurements of Leakage in Mechanical Face 
Seals . 

Presented at the Fourth Inti. Conference 
on Fluid Sealing held in Conjunction With 
24th ASLE Annual Meeting in Philadelphia, 
FICFS Preprint No. 19. 

32. Adams, M.L. and Colsher, R.J.(1969) 

An Analysis of Self, Energiged Hydrostatic 
Shaft Seals. 

Jour. Lub.Tech. ,Trans.ASME, Vol.91,No.4, 

pp. 658 . 

33. Nau, B.S. (1969) 

Film Cavitation observation in Face 
Seal . 

Presented at the Fourth International 
Conference on Fluid Sealing held in 
Conjunction with the 24th ASLE Annual 
Meeting in Philadelphia, FICFS, Preprint 
No.2 0. 


34. Lymer,A, (1969) 

An Engineering Approach to Selection and 
Application of Mechanical Seals. 

Presented at the Fourth International 
Conference on Fluid Sealing held in 
Conjunction with the 24th ASLE Annual 
Meeting, Philadelphia, FICFS Preprint No .2 5. 



151 


35. Orcutt.F.K. (1969) 

vin Investigation of the Operation and 
Failure of ^''■‘echanical Face Seals. 

Presented at the Fourth International 
Conference on Fluid Sealing, held in 
Conjunction with the 2 4th i'lSLE Annual 
Meeting, Philadelphia, FICFS Preprint No. 2 2 . 

36. Gibson, I.A Hooke, C.J. and O, Donoghue, J.P. (1970 ) 

Elas tohydrodynamic lubrication of 0-ring 
seals , 

Proc - Instn.Mech.Engrs . , Vol.l4, p.34 . 

37. Iny,E.H. (1971) 

The Design of Hydrodynamic Rubricated 
Seals With Predictable Operating Character- 
istics . 

5th Inti. Conf. on Fluid Sealing, Coventry, 
England, Paper HI. 

38. Iny,E.H. (1971) 

A Theory of Sealing With Radial Face Seals. 
Wear, Vol. 18,p.51. 

39. Iny, E.H., Stangham-Batch, B.A. and Winney, P.E.( 1-7 1 } 

Face Lubrication In Mechanical Face Seals. 
Inst.Mech.Engrs ., C 59/71, p.54. 

40. Field,G.J. and Nau, B.S ., (1972 ) 

An Experimental Study of Reciprocating 
R^Ibber Seals. 

Proc. Inst, of Mech.Engrs. Symposium on 
Elas tohydrodynamic Lubrication. 

<^■1. Lohou, J. and Godel,M. (1973 ) 

x'lngular Misalignment and Squeeze Film 
Effects In Radial Face Seals. 

6th Intl^' Conf. on Fluid Sealing, Munich, 
Published by BHRA, Cranfield, England, 

Paper B2 . 

42. Stangham-Batch, B.A. and Iny, S.H. (1973 ) 

A Hydrodynamic Theory of Radial Face Seals. 
Jour, of Mech.Engng.Sc. , Vol. 15, No. 1. 

43. Snapp,R.B. and Sasdelli, K.R. (1973 ) 

Performance Characteristics of a High 
Pressure Face Seal With Radially Converging 
Interface Shapes. 

Sixth Inti. Conf. on Fluid Sealing, Munich, 
Paper E4. 



152 


^ 4 , Metcalf e,R . ( 1973 ) 

Performance Analysis of Axisymmetric 
Flat Pace Mechanical Seals* 

Sixth Inti. Conf. on Fluid Sealing, 

Munich, Paper Dl. 

45. Mayer, E.(1973) 

Mechanical Seals . 

Iliffe Books Limited, London, American 
Elsevier, New York, II Edition, 250 pages. 

46. Zuk, J. (1974) 

Analytical Study of Pressure Balancing 
in Gas Film Seals. 

Trans .ASLE, Vol . 17,I.To.2 , p.97. 

47. Shapiro, W. and Colsher, R.{1974) 

Steady State and Dynamic Analysis of 
Jet Engine Gas -rubricated Shaft Seals. 
Trans. ASLE, Vol . 17, No .3 , p.l90. 

48. Haardt, R. (1974) 

Flow Consideration Aroxind Cavitation Area 
in Radial Face Seals. 

Proc. I Leeds-Lyon - Lyon Symp. on 
Tribology, Paper VIII ( i ) . 

^- 9 . Ruskell, L.E.C. and Westcotl,M, J. ( 1974 ) 

High Performance Reciprocating Seals for 
Aircraft Hydraulic System. 

Tribology, Vol. 5, p- 365. 

50. Haardt, R. and Godet, M. (1975) 

Axial Vibration of a Misaligned Face 
Seal Under a Constant Closure Force. 

Trans. ASLE, Vol.l8,No. l,p.55 • 

51. Griskin,E.N. (1975 ). 

The Effect of Dynamics on Fluid Flow 
in a Face Seal. 

Proc. 7th Inti. Conf. on Fluid Sealing, 
Paper B2 . 

52. Kupperman, D.S,(1975) 

Dynamic Tracking of Hon Contacting Pace 
Seals, 

Trans. ASLE, Vol.l8,No.4,p.306, 



153 


53. Field, G.J, and Mau, B.S.(l975) 

A Theoretical Study of the Blastohydrodynamic 
Lubrication of Reciprocating Rubber Seals. 

J . Lub . Tech, Trans .ASI IE, Vol * 18, p. 48 . 

54. Ludwig, L.P. (1976) 

Face Seal Lubrication I-Rroposed and 
Published Models* 

NASA, Technical Note i>iASA, Ti’JD-8101, 

55. Ludwig, L.P. and Allen G.P,(1976) 

Face Seal Liibrication II -Theory of 
Response to Angular Misalignment. 

NASA, Technical Note NASA, TND-8102 . 

56. Metcalf e,R. (1976) 

The use of Finite Element Deflection 
Analysis in Performance Predictions for 
End Face Seals. 

Presented at the Third Symposium 
on Engg. Application of Solid Mechanics, 
University of Ibronto. 

57. Li Chin-Hsiu (1976) 

Thermal Deformation in a Mechanical Face 
Seal . 

Trans .ASLE, Vol. 19, p, 146. 

58. Kilaparti, S.R. and Burton, R. A. (1976) 

Pressure Distribution for Patch-Like Contact 
in Seals with Frictional Heating, Thermal 
Expansion and Wear. 

J . Lub . Tech ., Trans . ASME, Vol. 98, No. 4, 
p.569. 

59. Ruskell,L.E.C, (1976 ) 

Reynolds Equation and Elastohydrodynamic 
Lubrication in Metal Seals. 

Proc. Royal Soc. London, A. 349, p»383 . 

6C , Lebeck,A .0 . , Teale J.L. and Pierce, R.E. (1977 ) 

Hydrodynamic Lubrication and Wear in 
Wavy Contacting Face Seals. 

Presented at the 1977 ASLE— ASME Lxibrication 
Conference, Kansa, City,ASME Paper No.77- 
Lub-18 . 



154 


61. Lebeck, A.O. (1977 ) 

A Study of Mixed Lubrication in 
Contacting Mechanical Face Seals. 
Presented at the Fourth Leeds -Lyon 
S^K^mposium on Lubrication, Lyon. 

62 . Lebeck, A . O , ( 1977 ) 

Mechanical Loading-Ti. Primary Source of 
of Waviness in Face Seals. 

Trans . ASLE, Vol .2 0, Mo . 3 , p . 195 . 

63. Lebeck, A.O., Teale, J.L. and Pierce, R.E, (1978 ) 

Hydrodynamic Lubrication with Wear and 
Asperity Contact in Mechanical Face Seals. 
Annual Report ME-36(78) OMR-414-1, 

Prepared for the Office of Naval Research 
under Contact No .ONR-N-00014-76-C-0071. 

64. Etsion, I. (1978) 

The Accuracy of Narrow Seal Approximation 
in Analysing Radial Face Seals. 

Trans .ASLE, Vol .2 2 , No .2 , p .2 08 . 

65. Etsion,!. (1978) 

Nonaxisymmetric Incompressible Hydrostatic 
Pressure Effects in Radial Face Seals. 

J . Lxib . Tech . , Trans .ASME, Vol . 100, No . 3 , 
p.379. 

66. Hughes, W.F., VTinowich, N.S,, Birchak, M.j. and 
Kennedy, W.C. (1978) 

Phase Change in Liquid Face Seals. 

J . Li±) . Tech • , Trans . ASl-iE , Vol . 100 , p . 7 4 . 

67. Patir, N, and Cheng, H.S. (197 8) 

An Average Flow Model for Determining 
the Effects of Three Dimensional Roughness 
on Partial Hydrodynamic Duijrication . 

J . Lub . Tech . , Trans .ASME, Vol . 100, p . 12 . 

68. Metcalfe,R,, Pothier,M.E. and Rod, B.H. (1978 ) 

Diametral Tilt and Leakage of End 
Face ^eals 'With Convergent Sealing 
Gap. 

Proc.Sth Intl-Conf. on Fluid Sealing, 

BHR\, 



155 


69. Kaneta, M., Fukahori, M, and Hirano, F. (1978) 

Dynamic Characteristics of Face 
Seals . 

Proc.Sth Inti. Conf. on Fluid 
Sealing. Paper A2 . 

'70. Lebeck, A.O. (1979) 

A Mixed Friction Hydrostatic 

Mechanical Face Seal Model With 
Thermal Rotation and Wear. 

Presented at the 34th ASLE Annual 
Meeting in St. Louis, ASLE Preprint 
N0.7 9-AM-4C-3 - 

71. Lebeck, A.O. (1979) 

A Mixed Friction Hydrostatic 

Face Seal Model With Phase 
Change . 

J. Lub- Tech. , Trans .ASME, Paper No. 
79-Lub-5,p.l. 

72. Sharoni, A. and Etsion, I. (1979) 

Performance of End Face Seals With 

Diametral Tilt and Coning — 

/ Hydrodynamic Effects. 

Trans. ASLE, Vol.24,No,l,p. 61. 

73. Etsion, I., (1979) 

Performance of End Face Seals 

With Diametral Tilt and Conting — 

Hydrostatic Effects. 

Trans.ASLE, Vol.23 ,No.3 , p.279 . 

74. Etsion, I., (1979) 

Hydrodynamic Effects in a 
Misaligned P-.adial Face Seal. 

J.Lub. Tech., Trans. ASME, 

Vol.lOl, p, 2 83. 

75. Etsion, I. and Burton, R.A. (1979) 

Observation of Self Excited 
Wobbles in Face Seals. 

J.Lub. Tech., Trans. ASME, Vol.lOl, 
No. 4, p.52 6. 



156 


76. Etsion, I. (1979) 

Squeeze Film Effects in Radial Face 
Seals . 

ASME Paper No.79-Lub~10. 

77. Etsion,I.and Sharoni.A, (1979 ) 

The Effects of Coning on Radial 
Forces in Misaligned Radial Face Seals. 
7.SME Paper No.79-Lub~17 . 


78. Etsion^I. (1979 ) 

Radial Forces in a Misaligned Radial 
Face Seal- 

J . Iiub - Tech. , Trans .ASME, Vol . lol/No . 1 , 

p.81. 

79. Ruskell,L.E.C. (1979) 

The Elastohydrodynamic Ferformance 
of Liow-Friction, Zero-Leakage Metal 
Seals . 

J .Mech.Engg.Sc. / Vol. 2 l,No.4^ p.2 75 . 


00, Lebeck,A.O . ( 1980 ) 

A Test Apparatus for Measuring the 
Effects of Waviness in Misaligned 
Face Seals. 

Presented at the ASLE Annual Meeting, 
lAnahiem, 


01. Young, L.?.. (1980 ) 

Experimental Evaluation of a Mixed 
Friction Hydrostatic Face Seal Model. 

Masters Thesis, The University of New 
Mexico, Albuquerque, New Mexico, 

82. Etsion, I. and Dan, Y. (1981), 

An Analysis of Mechanical Face Seal 
Vibrations. 

J . Lub . Tech . , Trans ,A3I4E, Vol - lo3 , p, 42 8 . 

03. Walowit, J .A . and Pinkus, 0,(1982 ) 

Iinalysis of Face Seals With Shrouded ockets. 
J. Lub. Tech, Trans. ASME, Vol.lo4,No.2 , 

p,2 62 . 



157 


84. Suyanami,T. Masuda, T., Oishi, N. and Shimazu.T. (1982 ) 

Study of Thermal Behaviour of Large Seal 
Ring. 

J.Lub.Tech. , Trans .ASME, Vol . 104 ^p* 449 . 

85 . Krauter^A . I . ( 1982 ) 

Measurement of Oil Film Thickness for 
Application to Elastomeric Stirling 
Engine Rod Seals. 

J.Lub. Tech- , Trans .ASME/ Vol . lo4/p. 455 • 

36. Etsion, I. (1982 ) 

Dynamic ivnalysis of Noncontacting Face Seals. 
J . Lub . Tech . , Trans - ASHE, Vol . lo4 / p . 4 60 . 

37. Mak,W,C. and Conway, H-D. (1977) 

The Lubrication of a Long, Porous, 

Flexible Journal Bearing. 

Jr. Lub. Tech, Trans. ASME, Vol.99,p,<!i49, 

33. Sinha, P- and Singh, C. (1982 ) 

Lubrication of a Cylinder on a Plane 
With a Non -Newtonian Fluid considering 
Cavitation . 

J . Lub . Tech . , Trans .ASME, Vol . lo4, p . 168 . 

09. Thomas, B.D. and John, W.H. (1956) 

i-dvances in Chemical Engineering. 

/academic Press Inc., 1st ed.,Vol.l,New York, 
p . lo3 . 

90. Dowson, D., and Taylor, C.M.(1966) 

Elastohydrostatic Lubrication of Circular 
Plate Thrust Bearing. 

ASME Paper No. 6 6 -Lub -S -15 . 

91. Bird,R.B. , Stewart, W.E. and Lightf oot, E . N . (I960) 

Transport Phenomena, Wiley, 

92. Shukal,J.B. and Prakash,J. (1969 ) 

The Rheostatic Thrust Bearing Using a 
Power Law Fluid as Lubricant. 

Japanese Jour, of Applied Phys,,Vol.8, 
p.l557. 

93. Pinkus,0., and StemliGht,B, (1961) 

Theory of Hydrodynamic Imbrication, 

McGraw Hill, New York. 



